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Abstract 

The  optimized  design  methodology  implemented  in  this  research  employs  a  design 
methodology,  namely  Response  Surface  Methodology  (RSM),  that  is  relatively  new  to  the 
electromagnetics  field.  The  design  methodology  used  a  full  factorial  designed  experiment, 
a  Hybrid  Finite  Element  Method  (HFEM)  analysis  code,  a  cubic  spline  empirical  model,  an 
Analysis  of  Variance  (ANOVA)  test,  and  a  simplex  optimization  code  to  successfully  mea¬ 
sure,  predict,  analyze,  and  optimize  average  antenna  gain  and  antenna  input  impedance  for 
a  cavity-backed  microstrip  antenna  with  dielectric  overlays.  Using  the  antenna  operating 
frequency,  the  relative  permittivities  of  two  dielectric  superstrate  layers,  and  the  thickness 
of  a  single  dielectric  superstrate  layer  as  experiment  variables,  the  full  factorial  designed 
experiment  used  a  HFEM  analysis  code  to  determine  average  antenna  gain  and  antenna 
input  impedance  for  five  levels  of  each  experiment  variable.  Based  on  HFEM  results  from 
the  full  factorial  designed  experiment,  a  cubic  spline  empirical  model  successfully  predicted 
average  antenna  gain  and  input  impedance  values  for  variable  combinations  within  the  ex¬ 
periment  variable  limits.  The  ANOVA  test  determined  the  impact  each  design  variable 
introduced  on  the  average  antenna  gain  and  indicated  the  areas  of  greatest  concern  during 
manufacturing  and  testing  phases.  Using  simplex  optimization,  an  optimized,  realizable 
antenna  design  with  minimum  average  antenna  gain  above  the  antenna  bore  sight  and  a 
target  input  impedance  of  50  ohms  was  located  within  two  minutes  that  agreed  within 
approximately  4.0  percent  to  HFEM  values  determined  from  design  validation  tests,  while 
measured  results  from  the  constructed  optimized  antenna  design  closely  replicated  the  pre¬ 
dicted  antenna  patterns  and  the  predicted  input  impedance  value.  The  design  methodology 
eliminated  the  need  of  numerous  and  costly  antenna  builds  by  applying  straight-forward 
RSM  techniques  and  saved  countless  man  hours.  The  significant  time  savings  found  us¬ 
ing  this  optimized  design  methodology  demonstrates  the  power  of  RSM  and  motivates  its 
application  to  other  electromagnetic  design  problems. 
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OPTIMIZED  DESIGN  METHODOLOGY  OF  CAVITY-BACKED 


MICROSTRIP  ANTENNAS  WITH  DIELECTRIC  OVERLAYS 

L  Introduction 

Proven  through  the  years,  antenna  design  generally  follows  two  principle  paths:  either 
original  research  or  existing  antenna  scaling.  In  general,  both  antenna  design  approaches 
require  in-depth  modeling  and  testing  that  use  established  numerical  methods  such  as 
Method  of  Moments  (MoM)  or  Finite  Element  Method  (FEM)  validation  prior  to  antenna 
construction.  Following  the  construction  of  a  physical  antenna,  additional  verification  of 
the  antenna’s  operational  characteristics  using  antenna  chamber  measurements  usually 
completes  the  overall  antenna  design  procedure.  This  process  demands  large  amounts 
of  time  and  computational  effort  as  well  as  expensive  prototype  builds.  A  need  exists 
in  the  antenna  design  field  for  a  new  design  methodology  that  maintains  antenna  design 
accuracy  while  providing  fast,  accurate  antenna  designs.  Response  Surface  Methodology 
(RSM)  and  designed  experiments  offer  a  solution.  This  research  applies  the  principles  of 
RSM  and  develops  an  optimized  antenna  design  methodology  that  uses  Hybrid  Finite  Ele¬ 
ment  Method  (HFEM)  predictions,  designed  experiments,  empirical  modeling.  Analysis  of 
Variance  (ANOVA)  tests,  and  simplex  optimization.  The  design  methodology  is  illustrated 
by  applying  the  principles  of  dielectric  loading  to  a  cavity-backed  microstrip  antenna  and 
determines  an  optimized  antenna  design  within  minutes  rather  than  days  without  using 
older,  time  consuming  antenna  design  approaches. 

lA  Problem  Statement 

This  research  is  tasked  to  develop  an  optimized  antenna  design  methodology  that  uses 
RSM  techniques  and  to  illustrate  its  application  on  a  cavity-backed  microstrip  antenna  by 
employing  the  principles  of  dielectric  loading.  The  general  antenna  configuration  includes: 

•  A  rectangular  cavity  with  constant  dimensions, 
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•  A  rectangular  radiating  aperture  with  constant  feed  pin  location  that  is  immersed  in 
a  dielectric  substrate  with  constant  relative  permittivity, 

•  Two  dielectric  superstrate  layers. 

The  optimized  design  methodology’s  goal  is  an  optimized  antenna  design  in  minutes  rather 
the  than  days  typically  required  using  traditional  antenna  design  approaches. 

1.2  Assumptions 

Assumptions  prior  to  beginning  research  progress  include: 

1.  The  ability  to  accurately  predict  antenna  pattern  shape  and  antenna  gain  using  a 
HFEM  code, 

2.  The  ability  to  control  antenna  pattern  shape,  antenna  gain,  and  antenna  input  im¬ 
pedance  using  dielectric  overlays, 

3.  The  ability  to  sufficiently  sample  antenna  gain  and  antenna  input  impedance  values 
using  designed  experiments, 

4.  The  ability  to  empirically  model  average  antenna  gain  and  antenna  input  impedance. 

1.3  Scope 

While  the  design  methodology  developed  here  could  be  applied  to  any  electromag¬ 
netic  design  problem,  for  the  purposes  of  illustration,  only  a  cavity-backed  microstrip 
antenna  design  problem  will  be  examined.  The  methodology  design  problem  accomplishes 
the  following  objectives: 

1.  Develops  an  optimized  design  methodology  based  on  HFEM  code  results  describing 
a  cavity-backed  microstrip  antenna  with  dielectric  overlays. 

2.  Creates  an  empirical  model  that  accurately  describes  the  antenna  gain  results  and 
input  impedance  in  a  region  of  interest. 

3.  Validates  the  empirical  model  within  the  region  of  interest. 


1-2 


4.  Applies  the  principles  of  ANOVA  testing  to  determine  the  effects  of  experiment  vari¬ 
ables  on  antenna  gain  thereby  indicating  the  antenna  gain’s  sensitivity  to  variation 
in  experiment  variable  levels. 

5.  Applies  an  optimization  routine  to  the  empirical  model  to  locate  and  define  an  opti¬ 
mized  antenna  configuration  in  the  desired  region  of  interest, 

6.  Validates  the  empirical  model  based  optimized  antenna  design  by  comparing  pre¬ 
dicted  model  results  with  HFEM  code  results  and  constructed  antenna  measurements 
of  identical  antenna  configurations. 

1.4  Resources 

Resources  necessary  for  the  research  accomplishment  include: 

1.  An  HFEM  code  that  specifically  accommodates  the  analysis  of  cavity-backed  mi¬ 
crostrip  antennas, 

2.  A  computer  resource  capable  of  performing  large  numbers  of  HFEM  code  predictions 
within  a  short  period  time, 

3.  A  compact  anechoic  range  facility  for  designed  antenna  measurements. 

1.5  Overview 

This  document  contains  a  review  of  the  recent  significant  work  performed  in  the  area 
of  dielectric  loading,  the  theory  behind  RSM  and  designed  experiments,  the  development 
of  an  optimized  design  methodology,  and  the  optimized  design  methodology  applied  to  a 
realizable  antenna  design  with  both  predicted  and  measured  antenna  results.  Chapter  II 
reviews  current  literature  available  in  the  area  of  dielectric  loading  effects  of  microstrip 
antennas,  antenna  analysis  techniques  using  MoM  and  HFEM  codes,  and  an  original  at¬ 
tempt  of  microstrip  antenna  optimization  using  an  adapted  form  of  RSM  and  designed 
experiments.  Chapter  III  provides  the  theory  behind  the  optimized  design  methodology 
developed  later  in  Chapter  IV.  Chapter  IV  develops  the  optimized  design  methodology 
using  a  specific  antenna  configuration  and  validates  the  design  methodology  using  stan¬ 
dard  designed  experiments  techniques  and  HFEM  results.  Chapter  V  applies  the  design 
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methodology  previously  developed  in  Chapter  IV  to  a  realizable  antenna  design  using  a 
common  substrate  dielectric  constant  and  designs  an  optimized  antenna.  Chapter  V  then 
validates  the  designed  antenna  using  two  approaches:  HFEM  code  validation  and  measured 
antenna  parameter  validation.  Chapter  VI  draws  conclusions  from  the  results  found  in  the 
previous  chapters  and  recommends  areas  where  this  design  methodology  is  applicable. 
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IL  Background 


2, 1  Introduction 

The  optimized  design  methodology  developed  in  Chapter  IV  and  applied  in  Chap¬ 
ter  V  uses  the  principles  of  dielectric  overlays  on  a  cavity-backed  microstrip  antenna  to 
control  antenna  gain  and  input  impedance  with  the  added  effect  of  controlling  the  an¬ 
tenna  radiation  pattern.  Chapters  IV  and  V  directly  apply  the  dielectric  overlay  effects 
described  in  Section  2.2  to  an  optimized  design  methodology  development  by  using  an 
analysis  technique  similar  to  those  described  in  Section  2.3.  Section  2.4  introduces  an 
initial  optimization  example  applying  an  adapted  form  of  Response  Surface  Methodology 
(RSM)  and  designed  experiments  to  microstrip  antenna  design.  In  a  similar  fashion,  the 
optimized  design  methodology  of  Chapter  IV  applies  a  form  of  RSM  to  a  cavity-backed 
microstrip  antenna  with  dielectric  overlays  making  Sections  2.2,  2.3,  and  2.4  directly  ap¬ 
plicable. 

2^2  Dielectric  Overlay  Effects 

According  to  the  effects  of  Snell’s  Law,  image  theory,  and  boundary  conditions, 
dielectric  overlays  strongly  influence  antenna  frequency,  gain,  input  impedance,  input  im¬ 
pedance  bandwidth,  and  frequency  bandwidth  [1-4,6,11,12,17,18,21].  For  example,  an 
optimized  selection  of  superstrate  permittivity  or  permeability  accompanied  by  an  opti¬ 
mized  selection  of  superstrate  thickness  positively  affects  antenna  gain  while  negatively 
affecting  other  antenna  properties  [2, 11,21]. 

2.2. 1  Resonant  Frequency.  Using  the  general  microstrip  antenna  design  displayed 
in  Figure  2.1,  Bahl  demonstrated  that  loading  a  microstrip  antenna  with  a  dielectric  su¬ 
perstrate  altered  the  antenna’s  resonant  frequency  [3].  The  absolute  resonant  frequency 
change  increased  with  operating  frequency,  superstrate  permittivity,  and  the  superstrate 
thickness  where  the  fractional  resonant  frequency  change  relative  to  the  unloaded  case  can 
be  calculated  using  Equation  (2.1)  for  a  fixed  superstrate  thickness  [3]. 

A/r  _  frjd  =0)-  frjd) 
fr  fr{d  0) 


(2.1) 
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Figure  2.1  Microstrip  Antenna  With  Dielectric  Superstrate  Geometry 

In  Equation  (2.1),  fr{d  =  0)  and  fr{d)  represent  the  unloaded  resonant  frequency  case 
and  the  loaded  resonant  frequency  case  respectively.  Equation  (2.2)  provides  a  first  order 
resonant  frequency  change  in  terms  of  effective  permittivity  values 
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where  and  eg  represent  the  effective  permittivity  without  superstrate  loading  and  the 
effective  permittivity  with  superstrate  loading  respectively  for  a  fixed  superstrate  thick¬ 
ness.  Bahl  observed  that  the  antenna’s  operating  frequency  decreased  as  the  superstrate’s 
permittivity  increased  and  as  the  superstrate’s  thickness  increased  [3]. 


2,2,2  Antenna  Gain.  Printed  circuit  antennas  suffer  from  the  disadvantage  of 
low  antenna  gain.  Using  Green’s  functions  and  contour  integration,  Alexopoulos  derived 
that  dielectrically  loading  an  embedded  Hertzian  dipole  with  a  single  dielectric  superstrate 
affects  the  resulting  antenna  gain  [2],  Specifically,  the  proper  choices  of  dielectric  constants 
and  thicknesses  for  the  substrate  and  superstrate  layers  increased  the  antenna  gain  and 
led  to  the  development  of  two  resonance  conditions.  Tables  2.1  and  2.2  list  the  details  for 
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Superstrate  thickness 

0.25  A 

Substrate  thickness 

0.50  A 

Dipole  location 

Substrate  center 

Permittivity 

^superstrate  ^  ^substrate 

Table  2.1  First  Resonance  Condition 


Superstrate  thickness 

0.25  A 

Substrate  thickness 

0.25  A 

Dipole  location 

Substrate/superstrate  interface 

Permeability 

f^superstrate  ^  (^substrate 

Table  2.2  Second  Resonance  Condition 


each  condition.  Application  of  the  resonance  conditions  led  to  the  development  of  a  simple 
transmission  line  model  for  quick  and  accurate  gain  calculations. 

Extending  the  ideas  developed  by  Alexopoulos,  Yang  introduced  a  transmission  line 
model  describing  an  embedded  Hertzian  dipole  antenna  covered  with  multiple  superstrate 
layers  [21].  Figure  2.3  illustrates  the  transmission  line  model  describing  the  antenna  de¬ 
sign  of  Figure  2.2.  From  a  ray  optics  perspective,  rays  emanating  from  the  antenna  and 
incident  upon  the  interface  of  different  media  will  bend  according  to  Snell’s  Law.  With 
a  transmission  line  analysis,  Yang  found  that  the  proper  arrangement  of  multiple  layer 
parameters  bent  transmitted  rays  in  a  prescribed  direction.  The  antenna  design  used  by 
Yang  displayed  in  Figure  2.2  assumes  an  infinite  ground  plane,  k  superstrate  layers,  and 
infinite  dielectric  layers  in  length  and  width.  Using  the  resonance  conditions  established 
by  Alexopoulos  for  the  single  layer  superstrate  design,  Yang  noticed  that  when  using  con¬ 
figurations  of  multiple  superstrate  layers,  the  alternation  of  layers  of  large  permittivity 
with  layers  of  large  permeability  yields  desired  high  antenna  gain.  Specifically,  in  the  case 
of  the  first  resonant  condition,  the  even  numbered  layers  have  large  permittivity  values, 
and  the  odd  numbered  layers  have  large  permeability  values.  In  the  case  of  the  second 
resonant  condition,  the  even  numbered  layers  have  large  permeability  values,  and  the  odd 
numbered  layers  have  large  permittivity  values.  When  resonant  conditions  are  satisfied 
and  layer  parameters  are  met,  the  high  gain  condition  described  by  the  transmission  line 
model  corresponds  to  the  maximum  voltage  transfer  from  the  incident  voltage  wave  Vin  to 
the  line  voltage  of  the  substrate  Vi,  at  a  prescribed  incident  angle  6  [21].  For  the  first  res- 


2-3 


Figure  2.2  Hertzian  Electric  Dipole  (HED)  Embedded  In  Substrate  With  Multiple  Su- 
perstrates 

Vta  V,  V,.,  V3  V,  Vj 

O - O - O  O - O - O -  , 

^  0 

I  ^  ^  j  ^  I  ^  ^  I  ^  ^  I  ^  ^ 

z,  .... 

o - o - o  o - o - o - 

Figure  2.3  Equivalent  Transmission  Line  Model 


onant  condition  at  a  prescribed  incident  angle  0,  the  maximum  voltage  transfer  for  2k-\-l 


layers  is  described  by 


-  Vin 


ZclZczZch  •  *  *  ^c(2A;+l) 

ZcoZc2ZcA  *  *  *  Zc{2k) 


while 


1^1  -  Vi, 


ZclZcsZc^  •  ’  •  ^c(2fc-l) 

Zc2Zc4  *  *  *  Zc(2k) 


(2.3) 


(2.4) 


describes  the  maximum  voltage  transfer  for  2k  layers.  In  Equations  (2.3)  and  (2.4),  Zd  is 
the  characteristic  impedance  of  layer  i  for  i  =  0, 1, 2,  •  •  • ,  A;.  Likewise,  in  the  case  of  the 
second  resonant  condition  at  a  prescribed  incident  angle  9,  the  maximum  voltage  transfer 
for  2A;  +  1  layers  is  described  by 


Vl  =  Vin 


Zc2ZcaZcQ  *  *  *  Zf,^2k) 
ZczZdsZc7  •  *  •  Zc(2k+l) 


(2.5) 


while 


Zc2ZciZcQ  •  •  •  Zc{2k) 

17  rr  ry  ry 

^c0^c3^c5  •  *  *  ^c(2fc-l) 


(2.6) 
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Figure  2.4  Microstrip  Antenna  With  Superstrata  Layer 

describes  the  maximum  voltage  transfer  for  2k  layers. 

Using  a  Method  of  Moments  (MoM)  approach  as  opposed  to  a  transmission  line 
model,  Shen  devised  another  antenna  design  technique  for  antennas  with  dielectric  super- 
strates  [17].  Consider  the  antenna  displayed  in  Figure  2.4.  Shen  discovered  that  using  a 
superstrata  layer  with  a  quarter  wavelength  thickness,  a  high  relative  permittivity  in 
layer  three,  and  adjusting  t2  led  to  high  antenna  gain.  Shen  observed  that  the  antenna 
gain  first  increased  with  the  distance  to  a  maximum  value,  Gmax^  at  a  resonant  dis¬ 
tance,  t2r*  For  every  ts,  a  corresponding  t2r  and  Gmax  exists  that  are  functions  of  the 
ta  plate  thickness.  Generally,  as  increases,  t2r  decreases,  and  Gmax  occurs  when  ts  is 
The  thickness  ti  is  not  of  immense  importance  since  the  thickness  of  layer  2  may  be 
adjusted  to  account  for  negative  effects  introduced  by  layer  one.  For  initial  considerations 
and  assuming  no  superstrate  loading,  the  initial  antenna  design  defines  the  thickness  of 
layer  one.  In  summary,  the  key  resonance  condition  is  the  electrical  thickness  of  ts,  and 
the  final  resonance  condition  can  be  tuned  by  adjusting  the  value  of  ^2- 
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2.2.3  Frequency  Bandwidth.  A  major  limitation  in  antenna  designs  implementing 
dielectric  loading  is  the  frequency  bandwidth.  According  to  Shen  and  Yang,  antenna 
frequency  bandwidth  decreases  as  antenna  gain  increases  [17,21].  In  other  words,  the 
frequency  bcuidwidth  continues  to  decrease  with  the  addition  of  gain  enhancing  superstrate 
layers.  Shen  defines  the  half  power  frequency  bandwidth  as 

A/o.5G  _  (/2  —  /i)  |g=o.5G  max  (2.7) 

/o  /o 

where  G  is  the  antenna  gain  and  /  is  the  frequency. 


2.2.4  Input  Impedance  and  Input  Impedance  Bandwidth.  A  second  limitation  to 
the  application  of  dielectric  superstrates  used  for  gain  enhancement  is  the  input  impedance. 
Soares  noticed  that  in  an  antenna  configuration  with  a  single  dielectric  superstrate,  the  an- 
terma’s  input  impedance  dramatically  decreased  as  the  anterma’s  superstrate  permittivity 
increased  [18].  The  only  way  to  maintain  an  adequate  input  impedance  level  is  to  reduce 
antenna  gain  when  implementing  dielectric  superstrates  for  gain  enhancement  [17]. 

In  comparison  with  the  frequency  bandwidth,  the  addition  of  superstrate  layers  re¬ 
duces  the  antenna’s  input  impedance  bandwidth  from  original  levels  without  dielectric 
loading  and  causes  the  input  impedance  bandwidth  to  become  even  more  restrictive  than 
the  frequency  bandwidth  [17, 18].  Shen  defines  the  input  impedance  bandwidth  as 

A/o.69fl,.  _  if 2  -  fl)  \R=0.69Rr  /2  S') 

fr  ~  fr 

where  R  is  the  input  impedance  at  0.69i2r,  Rr  is  the  resonant  input  impedance,  and  /  is 
the  frequency  at  specific  R  values.  When  matched  at  resonance,  the  impedance  bandwidth 
using  the  impedance  coefficient,  0.69,  has  a  return  loss  below  -10  dB.  The  consequences 
of  dramatic  input  impedance  level  changes  demand  that  input  impedance  requirements  be 
maintained  throughout  the  entire  antenna  design  process. 
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2.3  Cavity  Effects 


In  the  case  of  microstrip  antennas,  the  frequency  bandwidth  limits  their  usefulness. 
Thick  substrate  layers  circumvent  this  problem  but  create  the  additional  problem  of  surface 
wave  propagation.  In  the  case  of  cavity-backed  microstrip  antennas,  the  microstrip  antenna 
employs  a  metallic  cavity  to  enclose  the  antenna  and  suppress  surface  waves  [6].  In  addition 
to  substrate  mode  elimination,  the  metallic  cavity  serves  as  a  heat  sink  and  improves  heat 
dissipation  [6] .  Analysis  of  cavity-backed  aperture  antennas  performed  using  the  Method 
of  Moments  (MoM)  and  a  Hybrid  Finite  Element  Method  (HFEM)  are  discussed  in  the 
following  sections. 

2.3.1  Method  Of  Moments  Analysis  Of  Cavity-Backed  Antennas.  Due  largely 
to  the  development  of  high  speed  computers  with  enormous  computational  capabilities, 
complex  antenna  designs  can  be  solved,  modeled,  and  optimized  by  the  implementation 
of  MoM.  In  their  analyses,  Biebl  and  Lee  et  al.  demonstrated  the  principles  behind  MoM 
analysis  of  cavity-backed  aperture  antennas  [4,11].  Using  the  aperture  antenna  with  multi¬ 
ple  layers  in  Figure  2.5,  Biebl  implemented  MoM  to  analyze  the  effects  of  multiple  dielectric 
overlays  on  an  aperture  antenna.  Lee  on  the  other  hand  studied  an  antenna  with  a  single 
dielectric  overlay  using  MoM.  Both  Biebl  and  Lee  divided  the  antenna  into  two  regions, 
the  cavity  region  and  the  superstrate  region.  In  the  cavity  region,  dyadic  Green’s  functions 
applied  in  the  spatial  domain  corresponded  to  the  physical  design  and  shape  of  the  cavity 
itself  [11].  Similarly,  in  the  superstrate  region,  dyadic  Green’s  functions  applied  in  the 
spectral  domain  accounted  for  the  radiating  field  that  traveled  through  the  dielectric  over¬ 
lays  and  into  free  space  beyond.  By  modeling  the  superstrate  region  as  a  series  connection 
of  transmission  lines  and  the  cavity  region  as  a  waveguide,  MoM  solved  the  respective  re¬ 
gional  integral  equations  based  on  dyadic  Green’s  functions  and  determined  the  resulting 
solution  by  applying  a  weighted  series  summation.  In  the  superstrate  region,  the  dyadic 
Green’s  function  accounted  for  wave  impedance’s  for  both  the  transverse  magnetic  mode 
and  the  transverse  electric  mode,  and  in  the  cavity  region,  the  wave  propagation  constant, 
the  attenuation  factor,  and  ohmic  power  losses.  The  MoM  solution  was  only  valid  for 
operating  frequencies  greater  than  the  cavity  cutoff  frequency.  With  MoM,  the  double  in- 
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Figure  2.5  Geometry  Of  Cavity-Backed  Aperture  Antenna  With  Overlays 


tegral  dyadic  Green’s  function  was  solved  directly  without  the  introduction  of  asymptotic 
approximation  techniques. 

Using  the  MoM  technique,  the  effects  of  dielectric  loading  of  cavity-backed  antennas 
coincided  well  with  dielectric  overlay  effects  on  printed  circuit  antennas.  Lee  observed  that 
as  the  permittivity  of  the  single  superstrate  layer  increased,  the  resonant  frequency  and  the 
input  impedance  of  the  cavity-backed  aperture  antenna  decreased.  In  the  same  fashion,  as 
the  thickness  of  the  single  superstrate  layer  increased,  the  resonant  frequency  and  the  input 
impedance  of  the  cavity-backed  aperture  antenna  decreased.  The  superstrate  permittivity 
and  thickness  effects  of  the  cavity-backed  aperture  antenna  are  dual  effects  of  the  printed 
circuit  antenna  with  dielectric  overlays. 


2.3.2  Hybrid  Finite  Element  Method  Analysis  Approach.  Cheng  et  al.  [6]  applied 
a  hybrid  analysis  technique  to  a  complex  cavity-backed  microstrip  antenna  with  dielectric 
overlays  that  emphasized  the  positive  attributes  of  both  MoM  and  Finite  Element  Method 
(FEM)  principles  known  as  the  Hybrid  Finite  Element  Method  (HFEM)  and  proved  supe¬ 
rior  to  the  pure  MoM  analysis  procedure.  As  in  the  MoM  procedure,  the  hybrid  technique 
established  two  analysis  regions,  the  superstrate  region  and  the  cavity  region.  Although 
computationally  efficient  in  solving  spatially  unbounded  problems,  MoM  is  applicable  only 
to  a  few  specific  cavity  configurations  with  known  Green’s  functions.  Similarly,  FEM  is 
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useful  only  for  bounded  regions  of  known  volume  and  is  unreliable  for  unbounded  volumes. 
Consequently,  Cheng  applied  MoM  to  solve  the  dyadic  Green’s  functions  for  the  spectral 
domain  in  the  superstrate  region  and  FEM  to  solve  the  electromagnetic  field  distributions 
due  to  the  physical  complexities  of  the  cavity  region.  FEM  divided  the  cavity  configuration 
into  discrete  volumetric  portions  and  analyzed  each  individual  portion  separately.  Upon 
completion  of  the  individual  portion  analysis,  the  FEM  code  summed  the  individual  elec¬ 
tromagnetic  solutions  and  determined  the  net  cavity  solution.  The  overall  antenna  solution 
combined  both  regional  solutions  obtained  from  MoM  and  FEM  and  satisfied  the  tangen¬ 
tial  boundary  conditions.  Specifically,  the  tangential  field  components  resulting  from  the 
MoM  spectral  solution  of  the  superstrate  region  equaled  the  tangential  field  components 
resulting  from  the  FEM  spatial  solution  of  the  cavity  region  at  the  boundary  separating 
the  two  regions. 

2,4  Application  Of  Designed  Experiments,  Antenna  Modeling,  and  Parameter  Optimiza¬ 
tion 

Due  to  the  effects  and  limitations  of  dielectric  overlays,  the  design  of  microstrip 
antennas  with  single  or  multiple  dielectric  overlays  requires  numerous  tests,  builds,  and 
measurements,  A  new  design  approach  used  by  Moosbrugger  et  al,  [15]  implemented  a 
design  methodology  composed  of  several  prototype  antennas,  a  Global  Response  Surface 
(GRS)  model,  and  a  parameter  optimization  routine.  Prom  the  antenna  design  in  Fig¬ 
ure  2.6,  the  design  methodology  used  antenna  frequency,  resonant  patch  length  L^,  and 
dielectric  layer  separation  S  as  experimental  input  parameters  and  driven  patch  input 
impedance  as  the  experimental  measured  response.  The  remaining  physical  parameters 
labeled  in  Figure  2.6  remained  fixed  values  and  are  listed  in  Table  2.3.  Using  five  resonant 
patch  lengths,  four  layer  separation  distances,  and  31  frequencies,  the  design  methodology 
created  a  5  x  4  x  31  full  factorial  designed  experiment  where  full  factorial  designed  exper¬ 
iments  measure  the  system  response  at  each  possible  combination  of  the  input  variables. 
In  this  case,  the  full  factorial  designed  experiment  measured  the  antenna  input  impedance 
for  each  of  the  resulting  620  (5  x  4  x  31)  antenna  configurations.  Table  2.4  lists  the  ranges 
for  the  experimental  variables.  Based  upon  the  completed  impedance  measurements  set. 
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Top  View  of  Resonant  Patch 


Substrate,  £2 


Top  View  of  Driven  Patch 


Substrate, 


Air  or  Dielectric 


Ground  Plane 


Figure  2.6  Multi-layer  Microstrip  Patch  Antenna 
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Parameter 

Fixed  Value 

ei 

2.33 

^2 

2.55 

1.2mm 

h 

0.8mm 

Ld 

38mm 

Wd 

50mm 

Wr 

50mm 

Table  2.3  GRS  Model  Fixed  Parameters 


Parameter 

Symbol 

Range 

Resonant  patch  length 

Lf 

38.1mm  to  46.1mm 

Dielectric  layer  separation 

s 

5.1mm  to  10.1mm 

Frequency 

f 

2.ZGHZ  to  2.^GHz 

Table  2.4  GRS  Variable  Ranges 


the  implementation  of  a  GRS  model  used  a  cubic  spline  algorithm  to  predict  driven  patch 
impedance  values  at  various  variable  value  combinations  within  the  experimental  limita¬ 
tions  of  frequency,  layer  separation,  and  resonant  patch  length.  The  additional  application 
of  a  simplex  method  optimization  routine  allowed  the  GRS  model  to  locate  optimum  im¬ 
pedance  and  corresponding  variable  values  in  specific  frequency  ranges.  Table  2.5  describes 
an  optimized  antenna  configuration  completed  using  this  design  methodology  where  each 
variable  value  is  within  its  range  limits.  A  fabricated  antenna  constructed  from  the 
Table  2.5  design  agreed  extremely  well  with  predicted  GRS  impedance  and  return  loss 
values.  The  design  methodology  described  here  required  less  than  two  minutes  to  locate 
an  optimized  antenna  design  since  the  GRS  cubic  spline  antenna  model  required  less  than 
one  second  for  impedance  determination.  Moosbrugger  et  al.  demonstrated  that  under 
controlled  conditions  antenna  design  speed  is  greatly  improved  while  maintaining  design 
accuracy. 


Parameter 

Symbol 

Value 

Resonant  patch  length 

Ijf 

41.30mm 

Dielectric  layer  separation 

s 

7.12mm 

Frequency 

f 

2.5GHz  to  2.7GHz 

Driven  patch  input  impedance 

^opt 

20Q 

Table  2.5  Optimum  GRS  Antenna  Design 
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2.5  Summary 

The  effects  of  dielectric  superstrate  layers  on  printed  circuit  antennas  as  well  as 
cavity-backed  antennas  must  be  included  in  the  final  antenna  configuration.  The  selection 
of  dielectric  superstrate  and  substrate  permittivity  and  thickness  values  strongly  affect 
antenna  resonant  frequency,  gain,  input  impedance,  and  frequency  bandwidth  in  differ¬ 
ent  ways.  Proven  by  transmission  line  theory,  MoM  techniques,  and  HFEM  procedures, 
the  proper  selection  of  dielectric  constants  and  thickness  for  superstrate  layers  improves 
antenna  gain  while  decreasing  the  antenna’s  resonant  frequency,  input  impedance,  and  fre¬ 
quency  bandwidth.  To  reduce  the  difficulties  encountered  during  the  antenna  design  phase, 
a  new  antenna  design  methodology  was  introduced  that  corresponds  well  with  established 
antenna  design  methods  and  requires  only  seconds  for  completion.  The  design  methodology 
developed  in  Chapter  IV  incorporates  the  basic  design  methodology  outlined  in  Section  2.4 
but  uses  rigorous  HFEM  code  solutions  instead  of  multiple,  expensive  antenna  builds  to 
determine  antenna  response  data. 
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Ill,  Research  Methodology 


Response  Surface  methodology  (RSM)  comprises  a  set  of  statistical  techniques 
for  empirical  model  building  and  model  exploitations.  By  careful  design  and 
analysis  of  experiments,  it  seeks  to  relate  a  response,  or  output  variable,  to  the 
levels  of  a  number  of  predictors,  or  input  variables,  that  affect  it  [5]. 

In  the  development  of  an  optimized  design  methodology  for  cavity-backed  microstrip 
antennas  with  dielectric  overlays,  this  chapter  develops  the  theory  of  RSM  including  de¬ 
signed  experiments,  system  modeling,  and  system  optimization  used  throughout  this  re¬ 
search. 

3. 1  RSM  Tasks 

The  opening  quote  from  Box  and  Draper  outlines  the  key  RSM  design  methodology 
tasks  that  include: 

1.  Design  a  series  of  experiments  that  yields  reliable  measurements  of  the  system  re¬ 
sponse  of  interest  in  a  region  under  consideration. 

2.  Determine  an  empirical  model  that  accurately  describes  the  system  response. 

3.  Analyze  the  system  response  in  relation  to  the  effects  introduced  by  input  variables 
and  interactions  among  variables. 

4.  Search  for  the  optimal  levels  of  the  input  variables  that  produce  a  desired  system 
response. 

When  properly  completed,  the  four  tasks  produce  an  optimized  design  methodology  yield¬ 
ing  reliable  predicted  system  responses  over  the  region  of  interest. 

3.2  Experiment  Development 

In  response  to  the  first  RSM  task,  the  establishment  of  a  designed  experiment.  Box 
and  Draper  list  five  specific  initial  questions  that  include  [5]: 

1.  Which  variables  should  be  studied? 

2,  How  shall  qualitative  variables  be  chosen? 


3-1 


3.  At  which  levels  of  a  given  input  variable  should  experiments  be  run? 

4.  What  experimental  arrangement  should  be  used? 

5.  How  should  the  response  be  measured? 

This  section  addresses  each  question  and  provides  solutions  that  directly  apply  to  the 
development  of  an  optimized  design  methodology. 

3.2.1  Input  Variables.  The  first  question  of  Section  3.2  seeks  the  compilation  of 
a  set  of  applicable  input  variables  that  directly  impacts  the  system  under  consideration. 
The  first  attempt,  usually  a  brainstorming  session  with  all  interested  parties  followed  by 
a  variable  screening  procedure,  applies  working  knowledge  of  the  system  to  develop  a 
complete  variable  list  [5].  The  system  under  consideration  described  by 

r7  =  /(ei,6,---,a)  (3.1) 

contains  the  function,  /,  with  k  input  variables,  The  function  ,  /,  often  contains  a  large 
group  of  input  variables  with  strong  impacts  as  well  as  variables  with  weak  impacts.  The 
proper  selection  of  input  variables  found  here  helps  develop  later  an  empirical  system  model 
that  emphasizes  strong  input  variables  over  weaker  input  variables  thereby  reducing  the 
empirical  model  complexity  [7].  Throughout  the  experiment,  the  experimenter  must  remain 
flexible  and  continuously  review  the  variable  selection  in  the  event  previously  unconsidered 
variables  strongly  influence  system  responses. 

3.2.2  Variable  Screening  Procedure.  The  second  question  of  Section  3.2  compli¬ 
ments  the  selection  of  system  input  variables  by  providing  a  means  of  eliminating  nonessen¬ 
tial  variables  with  weak  system  impacts  thereby  keeping  the  number  of  useful  variables  to 
a  manageable  number.  As  in  the  variable  selection  stage,  the  first  attempt  of  list  reduction 
usually  applies  system  operating  knowledge.  Technical  knowledge  given  by  experts  often 
provides  the  simplest  approach  towards  reducing  the  variable  list  to  critical  factors,  and 
expert  technical  knowledge  potentially  uncovers  variables  overlooked  during  brainstorming 
sessions.  The  second  variable  list  reduction  technique,  a  two-level  full  factorial  designed 
experiment,  tests  the  effects  of  input  levels  and  their  functional  impacts  and  quickly  decides 


3-2 


the  necessity  of  retaining  particular  variables  in  the  experiment  [5].  Discussed  in  a  later 
section,  full  factorial  designed  experiments  provide  an  excellent  experimental  design  for 
the  gathering  of  system  responses  in  a  systematic  fashion.  Here,  in  the  variable  screening 
process  the  two-level  full  factorial  designed  experiment  tests  the  system  response  at  two 
separate  levels  of  each  variable  used  in  the  experiment.  For  a  total  of  n  variables  in  the  ex¬ 
periment,  all  possible  variable  combinations  are  implemented  and  their  responses  recorded 
for  a  total  of  2"  possible  variable  combinations  and  measured  responses.  To  shrink  the 
size  of  the  screening  experiment,  a  two-level  full  factorial  designed  experiment  of  only  the 
variables  in  doubt  could  be  performed  to  discover  variable  relationships  to  experimental 
system  responses. 

3.2.S  Variable  Levels.  The  third  question  of  Section  3.2  seeks  an  appropriate 
number  of  variable  levels  and  preludes  the  experimental  design  selection  used  for  gathering 
system  response  data.  In  addition,  the  number  of  variable  levels  chosen  here  directly 
impacts  both  the  empirical  model  accuracy  used  to  predict  system  responses  and  the  time 
required  to  gather  the  entire  set  of  measured  system  responses.  In  particular,  as  the  number 
of  variable  levels  increase,  the  size  of  the  experiment  increases,  and  the  time  necessary 
for  experiment  completion  drastically  increases.  On  the  other  hand,  the  accuracy  of  an 
empirical  system  model  based  on  an  experiment  including  three  or  more  variable  levels 
per  variable  far  exceeds  the  accuracy  of  a  similar  empirical  model  that  describes  the  same 
system  and  uses  only  two  variable  levels  [14].  The  mathematical  model  using  a  greater 
number  of  variable  levels  precisely  determines  system  responses  due  to  variable  changes 
more  accurately  than  models  using  fewer  variable  levels.  An  experimenter  must  reach  a 
decision  concerning  the  number  of  variable  levels  and  the  time  expense  necessary  for  an 
accurate  experimental  result. 

3.2.4  Experimental  Designs.  In  response  to  the  fourth  question  of  Section  3.2, 
two  common  experimental  designs  used  in  RSM  experiments  include  the  full  factorial 
designed  experimental  and  the  central  composite  designed  experiment.  Of  these  designs, 
the  central  composite  design  remains  the  most  widely  used  despite  numerous  advantages 
offered  by  the  full  factorial  design  [14]. 
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Variable  1 

Variable  2 

Response 

la 

2a 

1 

la 

2b 

2 

la 

2c 

3 

lb 

2a 

4 

lb 

2b 

5 

lb 

2c 

6 

Ic 

2a 

7 

Ic 

2b 

8 

Ic 

2c 

9 

Table  3.1  Two  Factor,  Three  Level,  Full  Factorial  Design 

3.2,4- 1  Factorial  Designed  Experiments.  Full  factorial  designed  exper¬ 

iments  are  detail  oriented  and  highly  accurate  regardless  of  system  response  smoothness 
or  roughness.  Furthermore,  full  factorial  designs  consider  the  system  response  for  each 
input  variable  level  as  well  as  responses  caused  by  variable  interaction  by  investigating 
system  responses  for  all  possible  variable  level  combinations.  For  example,  in  the  case 
of  the  system  described  in  Table.  3.1  with  two  input  variables  using  three  levels  each,  a 
complete  full  factorial  experiment  compiles  a  total  of  3^  meajsured  system  responses,  one 
for  each  variable  combination,  full  factorial  designs  are  extendible  to  any  number  of  input 
variables  with  the  unfortunate  cost  of  an  increased  number  of  measured  responses  and  an 
increased  time  requirement  for  experiment  completion  [13].  Montgomery  lists  several  ad¬ 
vantages  of  full  factorial  designed  experiments  over  all  other  forms  of  experimental  designs 
that  include  [14]: 

1.  More  efficient  than  a  one  factor  at  a  time  experiment, 

2.  Necessary  when  variable  interactions  are  present  to  avoid  misleading  conclusions, 

3.  Allows  the  effects  of  a  factor  to  be  estimated  at  several  levels  of  the  other  factors, 
yielding  conclusions  valid  over  a  range  of  conditions, 

4.  Allows  multitudes  of  comparisons  to  be  performed  and  so  facilitate  model  creation 
and  criticism, 

5.  Provides  highly  efficient  estimates  of  constants  whose  variances  are  small, 

6.  Leads  to  simple  calculations. 
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Variable 

One 

Variable 

Two 

la 

lb 

Ic 

2a 

lai2ai 

Ia22a2 

lbi2ai 

Ib22a2 

lci2ai 

Ic22a2 

\ajn2CLjn 

^bjri2cbjn 

lCm2n7n 

2b 

lai26i 

102262 

lbi2bi 

162262 

lci26i 

1C2262 

15m  25m 

lCm*^bjn 

2c 

lai2ci 

Ia22c2 

15i2ci 

1522c2 

lci2ci 

1C22c2 

^bfYi2Cfn 

lCyyj,2C77i 

Table  3.2  Pull  Factorial  Design  with  Two  Variables,  Three  Levels,  and  m  Replicates 

In  general,  m  replicates  of  the  experiment  are  performed  for  nondeterministic  systems 
providing  m  measured  system  responses  for  each  variable  combination  [14].  Assuming 
system  responses  with  normally  distributed  random  systematic  error,  m  replicates  provide 
a  means  of  determining  the  expected  system  response  for  each  variable  level  combination. 
Table  3.2  indicates  a  typical  full  factorial  design  using  m  replicates  with  two  variables  each 
having  three  levels. 

Central  Composite  Designed  Experiments.  Central  composite  de¬ 
signed  experiments  establish  strategic  measurement  locations  by  enclosing  the  system  re¬ 
sponse  space  with  system  measurements  and  by  including  repeated  system  measurements 
central  to  the  system  space.  Figure  3.1  represents  a  typical  central  composite  design  used 
to  gather  system  response  measurements  using  only  three  input  variables.  Central  compos¬ 
ite  designs  strongly  reduce  the  required  number  of  measurements  for  systems  with  smooth 
responses  and  are  generally  used  to  create  second  order  empirical  models  [13].  Therefore, 
central  composite  designs  offer  little  help  in  modeling  systems  with  sharp  transitions  or 
discontinuities  since  an  insufficient  number  of  measurements  are  available  to  accurately 
characterize  the  system  performance  throughout  the  entire  range  of  each  input  variable. 
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Variable  2 


0 


Figure  3.1  Central  Composite  Design  For  k  =  3 

Mason  lists  the  steps  needed  to  construct  a  central  composite  designed  experiment  with  k 
input  variables  [13]. 

1.  Construct  a  complete  2*'  full  factorial  designed  experiment  with  high  and  low  variable 
levels  depending  on  the  need  for  efficiency  and  the  ability  to  ignore  system  variable 
interaction. 

2.  Add  2k  axial  test  points  along  the  coordinate  axes. 

3.  Add  multiple  repeat  observations  at  the  experimental  design  center. 

3.2.5  System  Response  and  Measurement.  The  fifth  qnestion  of  Section  3.2 
determines  the  desired  system  response  and  the  means  by  which  the  response  is  measured. 
The  system  response,  often  called  the  main  effect,  is  generally  the  easiest  to  determine  since 
the  main  effect  is  usually  the  reason  for  conducting  the  experiment.  Once  again,  system  or 
functional  knowledge  is  useful  in  the  determination  of  an  appropriate  main  effect.  In  the 
case  of  antenna  design,  antenna  gain  or  input  impedance  are  excellent  main  effect  choices. 
Furthermore,  in  the  case  of  a  cavity-backed  microstrip  antenna  with  dielectric  overlays, 
the  Hybrid  Finite  Element  Method  (HFEM)  is  the  method  of  choice  to  predict  antenna 
responses.  The  HFEM  method  is  an  adapted  version  of  the  Finite  Element  Method  (FEM) 
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that  applies  boundary  conditions  linking  internal  antenna  cavity  space  to  external  cavity 
space. 


3.2.5. 1  Hybrid  Finite  Element  Method,  According  to  Jin,  the  HFEM 
method  incorporates  the  advantages  of  integral  equation  methods  and  FEM  methods 
to  solve  electromagnetic  problems  composed  of  bounded  and  unbounded  regions  such  as 
cavity-backed  aperture  antennas  [10].  The  general  principle  of  the  technique  introduces  a 
fictitious  boundary  separating  the  cavity  region  from  the  unbounded  region.  Interior  to 
the  cavity,  the  HFEM  method  formulates  electric  and  magnetic  field  solutions,  whereas  in 
the  exterior  region  HFEM  uses  integral  equations  to  determine  field  solutions.  At  the  ficti¬ 
tious  boundary,  field  continuity  conditions  connect  the  fields  in  the  two  regions,  leading  to 
a  coupled  system  for  the  solution  of  interior  and  exterior  fields.  In  the  external  region,  the 
HFEM  method  applies  the  well  known  vector  wave  equation,  the  dyadic  Green’s  function, 
and  the  Sommerfield  radiation  condition  to  determine  electric  and  magnetic  field  quantities 
at  a  distance  r  from  the  fictitious  boundary  as  well  as  located  on  the  boundary.  Since  the 
tangential  electric  and  magnetic  fields  must  be  continuous  across  the  boundary,  field  quan¬ 
tities  just  above  the  cavity  boundary  equal  those  just  below  the  boundary.  Section  3.2. 5. 2 
develops  the  general  FEM  approach  used  within  the  cavity. 

3. 2. 5. 2  Finite  Element  Method,  This  section  follows  the  general  FEM  de¬ 
velopment  used  by  Jin  [10].  FEM  analysis  replaces  an  entire  spatial  domain  by  a  number  of 
sub  domains  in  which  an  unknown  function  is  represented  by  simple  interpolation  fxmctions 
with  unknown  coefficients.  FEM  analysis  includes  the  four  basic  steps: 

1.  Discretization  of  the  domain, 

2.  Selection  of  interpolation  functions, 

3.  Formulation  of  system  equations, 

4.  Solution  of  the  system  of  equations. 

The  first  and  perhaps  the  most  important  step,  discretization  of  the  domain  fl,  affects 
computer  storage  requirements,  computational  time,  and  accuracy  of  the  numerical  re¬ 
sults.  Depending  on  the  domain  shape,  FEM  partitions  the  domain  using  common  ele- 


3-7 


ment  shapes  such  as  line  segments,  triangles,  triangular  prisms,  or  rectangular  bricks  and 
seeks  to  find  individual  solutions  relative  to  each  discrete  element.  The  second  step  se¬ 
lects  an  interpolation  function  that  provides  an  approximation  of  the  unknown  solution 
within  each  element.  The  interpolation  functions  of  each  element  are  usually  selected  to 
be  a  polynomial  of  first  or  second  order  since  higher  order  polynomials  usually  result  in 
a  more  complicated  composite  interpolation  function  formulation.  The  third  step  formu¬ 
lates  a  system  of  equations  by  using  the  Rayleigh-Ritz  or  the  Galerkin  methods  and  then 
sums  the  elemental  equations  over  all  elements  to  obtain  the  final  form  of  the  system  of 
equations.  Next,  Neumann  and  Dirichlet  boundary  conditions  are  applied  to  each  discrete 
element  thus  completing  the  system  of  equations.  The  final  step,  solving  the  system  of 
equations,  provides  an  approximation  to  the  fields  in  the  problem  domain  ft  for  which  the 
desired  antenna  parameters  such  as  gain,  input  impedance,  and  radiation  patterns  can  be 
measured. 

5.5  System  Modeling 

The  second  task  in  the  design  methodology  development,  system  response  modeling, 
allows  the  experimenter  to  gain  an  understanding  of  the  system  under  study  and  provides 
a  means  to  predict  system  responses  for  differing  input  variable  levels.  A  model  of  the 
system  response,  77,  composed  of  $  input  variables  and  0  constants  takes  on  the  general 
form 

T]  =  -\-e  (3.2) 

^  =  [^1)^2,  •  •  •  >Cfe]  (3-3) 

©  =  (3.4) 

where  the  error,  e,  is  desired  to  be  as  small  as  possible  [5].  Mathematical  models  are  of 
two  forms,  the  mechanistic  and  the  empirical  models, 

3.3.1  Mechanistic  Modeling.  Box  and  Draper  refer  to  the  mechanistic  model  as 
the  true  system  functional  form  based  on  initial- value  and  boundary- value  problems  such 
as  Neumann  and  Dirchlet  equations  [5] .  Mechanistic  models  use  known  physical  properties 


3-8 


to  scientifically  describe  the  complex  performance  and  responses  produced  by  the  system 
under  study.  Mechanistic  models  are  of  the  general  form 

V  =  (3.5) 

where  the  function  /  incorporates  all  known  physical  properties  and  provides  an  exact 
solution.  HFEM  is  an  example  of  a  mechanistic  model, 

3.3.2  Empirical  Modeling,  In  the  absence  of  scientific  or  physical  knowledge 
of  a  system  or  burdened  with  extremely  complex  mechanistic  models,  Box  and  Draper 
prefer  empirical  system  models  [5].  An  empirical  model  closely  approximates  the  true, 
mechanistic  model  and  takes  on  the  following  form: 

7?  =  /(#,©)  fa  p($,0)  +  e  (3.6) 

where  the  error,  e,  is  kept  as  small  as  possible.  Two  common  techniques  for  empirical 
model  building  include  least  squares  and  cubic  splines.  Least  squares  curve  fitting  is 
widely  understood  and  therefore,  is  not  developed  in  this  chapter.  The  interested  reader 
is  referred  to  experimental  design  and  linear  algebra  textbooks  such  as  Montgomery  and 
Strang  which  thoroughly  develop  least  squares  curve  fitting  [14, 19]. 

3.3.3  Cubic  Splines.  This  section  follows  the  cubic  spline  development  of  Gerald 

and  Wheatley  [9].  Cubic  splines  fit  functions  with  local  irregularities  without  violent 
misbehavior  and  retain  smoothness  in  regions  of  functional  smoothness.  In  the  development 
of  cubic  splines,  a  set  of  cubics  pass  through  a  set  of  measured  data  points  fitting  a  new 
cubic  to  each  interval  between  data  points.  For  the  pair  of  cubics  that  join  at  each  data 
point,  a  requirement  exists  that  the  slope  and  curvature  of  each  cubic  equal  the  slope  and 
curvature  of  the  joining  cubic  at  the  data  point.  The  cubic  for  the  ith  interval,  which  lies 
between  the  points  {xi,yi)  and  T/i+i)  takes  the  form 

y  ^  ai(x  -  Xi)^  +  bi{x  -  Xi)^  +  Ci{x  ~  Xi)  +  di  .  (3.7) 
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Since  it  fits  at  the  two  endpoints  of  the  interval, 


yi  —  ^ “f"  df 


(3.8) 


2/i+l  —  ^z)  “1”  ^i)  d"  ^(^i+1  ®i)  dj 

=  ajhf  +  bihi  +  Cihi  +  di 


(3.9) 


where 


hi  —  Xi^\  -  Xi  .  (3.10) 

To  fulfill  the  slope  and  curvature  requirement  for  joining  cubics,  Equations  (3.11)  and 
(3.12)  provide  the  first  and  second  derivatives  of  Equation  (3.7)  in  the  ith  interval. 


3ai{x  —  Xi)^  -1-  2bi{x  —  Xi)  +  Ci 

(3.11) 

y"  =  6ai(x  —  Xi)  -H  2bi 

(3.12) 

Writing  the  equations  in  terms  of  the  second  derivatives  of  the  interpolating  cubics  and 
letting  Si  and  5i+i  represent  the  second  derivatives  at  the  points  {xi,yi)  and  (xf+i,  j/j-i-i) 
simplifies  the  mathematical  procedure. 


Si  —  QOii^Xi  Xi^  21)^ 

=  2bi 


(3.13) 


5i+i  =  Qai{xi+i  -  Xi)  +  2bi 
=  6aihi  +  2bi 


(3.14) 


Prom  Equations  (3.13)  and  (3.14), 


bi  = 


(hi 


2  ’ 

~  Sj 

Ghi 


(3.15) 

(3.16) 
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Substituting  the  relations  for  Oj,  h,  and  di  from  Equations  (3.16),  (3.15),  and  (3.8)  into 
Equation  (3.9)  and  then  solving  for  Cii 


yi+i  =  ^  ^  +  2  ’ 

(3.17) 

Vi  *2hiSi  “h  hiSi-^\ 
hi  6 

(3.18) 

Invoking  the  equality  of  slopes  condition  for  two  cubics  that  join  at  {xi, 

T/i),  Equation  (3.11) 

becomes 

(3.19) 

In  the  previous  interval  from  Xi-i  to  a:*,  the  slope  at  {xi,yi)  becomes 

“  3uj — 1  («r j  Xi — i)  “1“  25j_i(irj  x^ — i)  d"  — x 

=  3ai—ihf_i  +  2bi—ihi—i  +  Cj—i  .  (3.20) 

Equating  Equations  (3.19)  and  (3.20)  and  substituting  for  a,  b,  c,  and  d  in  terms  of  S  and 
y  yields 

/  yi+l  yi  2h%Si  +  hiSi-\-l 
Vi  =  - 

=  3 

6 

Simplifying  Equation  (3.21)  produces 

hi-iSi-i  +  (2/ii_i  +  2hi)Si  +  hiSi+i  =  6  •  (3-22) 

Equation  (3.22)  applies  at  each  cubic  location  beginning  in  the  second  cubic  interval  and 
ending  in  the  next  to  last  cubic  interval  and  gives  n  —  2  equations  relating  the  n  values  of 
Si. 


2hi — xSi — 1  “h  hi — \Si 


(3.21) 
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Specifying  the  end  conditions  for  the  whole  curve  provides  the  two  equations  that 
involve  Si  and  5n.  Three  choices  exist  for  end  conditions  where  the  specification  of  which 
end  condition  is  somewhat  arbitrary.  The  first  end  condition  in  Equation  (3.23)  assumes 
that  end  cubics  approach  linearity  at  their  extremes. 

Si^Sn  =  0  (3.23) 


The  second  end  condition  in  Equation  (3.24)  assumes  that  end  cubics  approach  parabolas 
at  their  extremities. 

(5i  =  ^2),  (Sn  =  Sn-l)  (3.24) 


The  third  end  condition  in  Equations  (3.25)  and  (3.26)  treats  as  a  linear  extrapolation 
from  52  and  Ss  and  treats  Sn  as  a  linear  extrapolation  from  5n-i  and  5^-2 •  At  the  left 


end 


S2-S1  Ss-S2  {hi  +  h2)S2  ~  hiSs 

hi  h2  '  ^  h2 


At  the  right  end 


Sn  Sn—1  Sn—l  Sn—2 


h 


n—l 


hn—2 


5n- 


{hn—2  H“  hn—l)Sn—l  hn—lSn—2 

hn-2 


(3.26) 


The  first  condition,  often  called  a  natural  spline^  sometimes  flattens  the  curve  too  much  at 
the  ends.  The  third  end  condition  sometimes  gives  too  much  curvature  to  the  ends.  The 
best  curve  fit  is  given  when  one  knows  the  slope  of  the  function  at  the  endpoints. 

Writing  the  equations  for  52,  53,  . . .,  5^-1  from  Equation  (3.22)  in  matrix  form 
produces  Equation  (3.27).  In  Equation  (3.27),  n  —  2  equations  exist  with  n  unknowns  in 
the  5  vector.  The  selection  of  end  conditions  eliminates  two  unknowns,  5i  and  5^,  thereby 
reducing  the  5  vector  to  n  —  2  unknowns.  The  vector  5  may  now  be  solved.  After  finding 
values  for  5  in  each  interval,  values  in  the  ith  interval  for  Uj,  6^,  c^,  and  di  are  found  using 
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Equations  (3.16),  (3.15),  (3.18),  and  (3.8)  respectively. 


hi  2(^1  +  /i2)  ^2 

/i2  2(/i2  +  h%)  hz 


hn-2  2(/i„_2  +  hn-l)  /in-1 


-Si 

S2 

Ss 

Sn-l 

5„ 


2/3 -V2  _  V2-V1 

h2  h\ 

VA-yz  _  yz~y2 
hs  /l2 


j/n  S/n  — X  yn— 1  j/n— 2 

/in— 1  /in— 2 


(3.27) 


5.5..^  Model  Adequacy,  In  fitting  empirical  models,  model  adequacy  checks  re¬ 
quire  an  analysis  of  the  residuals  occurring  from  the  comparison  of  measured  and  modeled 
data.  The  major  plots  that  graphically  represent  the  fit  of  the  empirical  model  include  the 
normal  probability  plot  of  residuals,  the  plot  of  residuals  versus  fitted  values,  and  the  plot 
of  residuals  versus  individual  variable  parameters  [14].  Residuals  account  for  the  difference 
between  modeled  values  and  true,  expected  values. 

3.34^1  Normal  Probability  Plot  of  Residuals.  A  normal  probability  plot  of 
residuals  is  just  a  graph  of  the  cumulative  distribution  of  the  residuals  on  normal  proba- 
bility  paper,  that  is,  graph  paper  with  the  ordinate  scaled  so  that  the  cumulative  normal 
distribution  plots  as  a  straight  line  [14].  To  construct  a  normal  probability  plot,  arrange 
the  residuals  in  increasing  order  and  plot  the  fcth  probability  of  these  ordered  residuals 
versus  the  cumulative  probability  point  [14].  If  the  underlying  error  distribution  is  nor¬ 
mally  distributed  and  the  empirical  model  closely  fits  the  system  response  data,  the  normal 
probability  plot  will  resemble  a  straight  line  similar  to  the  plot  of  Figure  3.2  [14]. 

3. 3. 4.. 2  Plot  of  Residuals  versus  Fitted  Values.  A  plot  of  residuals  versus 
fitted  values  plots  the  empirically  modeled  values  versus  their  residual  values  on  a  standard 


Normal  Probability  Plot 


Figure  3.2  Normal  Probability  Plot  of  Least  Squares  System  Response  Model 

xy  graph.  For  an  empirical  model  containing  normally  distributed  random  systematic 
error,  residual  versus  fitted  value  plots  should  be  structure-less  thereby  indicating  that  the 
empirical  model  does  not  follow  any  specific  prediction  trend  [14].  For  example,  Figure  3.3 
represents  a  typical  residual  versus  fitted  value  plot. 

3, 3, 4-3  Plot  of  Residuals  versus  Individual  Variable  Parameters,  A  plot 
of  residuals  versus  individual  input  variables  plots  the  individual  variable  levels  versus 
empirical  model  residuals  thereby  graphically  indicating  the  variability  of  the  measured 
system  response  as  a  function  of  the  individual  input  variable  [14] .  The  amount  of  variation 
in  model  residuals  relating  to  specific  variable  levels  indicates  the  impact  that  each  variable 
level  has  on  the  prediction  accuracy  of  the  empirical  model.  Figure  3.4  represents  a  typical 
residual  versus  variable  parameter  plot  where  the  variable  level  near  3.45  introduces  less 
variation  in  the  empirical  model  accuracy  than  the  levels  near  2.85  and  3.75. 


3.4  Analysis  of  Variance 

The  third  task  in  the  design  methodology  development,  the  analysis  of  variable  effects 
on  the  measured  system  response,  uses  an  Analysis  of  Variance  (ANOVA)  test.  ANOVA 
tests  apply  to  both  full  factorial  designed  experiments  and  central  composite  designed 
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Variable 

One 

Variable 

Two 

1 

2 

b 

1 

yiii,  2/112 

■  ■  ■  )  2/lln 

2/121,  2/122  • 

•  •  •  ,  2/12n 

•  2/161 ,  2/162 

■  ■  ■  ,  2/l6n 

2 

2/211,  2/212 

•  •  •  ,  2/21n 

2/221,  2/222  • 

•••,  2/22n 

•  •  2/261 ,  2/262 

■  ■  ■  ,  2/26n. 

a 

Vail,  2/o12 

■  ■  ■  ,  Vain 

2/o21,  2/o22  • 

'  ■  ■  ,  2/a2n 

••  2/o61,  2/a62 

■  ■  '  ,  2/a6n 

Table  3.3  General  Two  Variable  Full  Factorial  Design 


experiments  with  any  number  of  experiment  input  variables  and  levels  with  the  single  re¬ 
quirement  that  at  least  two  replicates,  repeated  measurements,  are  taken  at  each  measure¬ 
ment  location.  This  section  closely  follows  Montgomery’s  development  of  the  two-factor 
ANOVA  test  that  describes  a  two- factor  full  factorial  designed  experiment  [14].  Using  a 
sum  of  squares  approach  and  a  Chi  Square,  distribution  test  statistic,  the  F  statis¬ 
tic,  ANOVA  tests  analyze  the  effects  introduced  by  experiment  variables  and  the  effects 
introduced  by  variable  interactions  on  the  measured  system  response. 

Using  the  general  case  of  the  two-factor  full  factorial  designed  experiment  with  a 
total  levels  of  variable  A,  b  total  levels  of  variable  J5,  and  n  replicates,  the  measured 
system  response,  Vijk^  for  the  ith  level  of  A,  the  jth  level  of  5,  and  the  kth  replicate 
appears  in  the  general  form  of  Table  3.3.  Equation  (3.28)  describes  the  measured  system 
response. 

Vijk  ~  Ti  +  Pj  +  (T/3)ij  -h  Cijk  (3.28) 

where: 

•  fjL  denotes  the  overall  mean  effect, 

•  Ti  denotes  the  effect  of  the  ith  level  of  variable  A, 

•  Pj  denotes  the  effect  of  the  jth  level  of  variable  5, 

•  {TP)ij  denotes  the  effect  of  interaction  between  and  /3j, 

•  Cijk  denotes  a  random  error  component. 
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34-1  Test  Hypotheses.  In  the  two-factor  ANOVA  test,  row,  column,  and  inter¬ 
action  effects  are  of  equal  interest.  Two  hypotheses  for  row  effects  state  that  either  the 
effects  for  all  row  levels  equal  zero,  or  at  least  one  row  effect  does  not  equal  zero. 

Hq:  Ti  =  T2  =  ---  =  Ta  =  0  (3.29) 

Hi  :  at  least  one  r*  0 

The  two  hypotheses  for  colmnn  effects  state  that  either  the  effects  for  all  column  levels 
equal  zero,  or  at  least  one  column  effect  does  not  equal  zero. 

Ho:  pi  =  P2  =  ---=l3b  =  0  (3.30) 

Hi  :  at  least  one  /  0 

Concerning  variable  interaction,  two  hypotheses  state  that  either  the  effects  for  all  row 
column  interaction  levels  equal  zero,  or  at  least  one  row  column  interaction  effect  does  not 
equal  zero. 

^0  :  =  0  h  3  (3.31) 

Hi  :  at  least  one  {rP)ij  ^  0 

The  hypotheses  seek  to  determine  whether  experiment  variables  influence  the  system  re¬ 
sponse  or  should  be  excluded  from  the  entire  experiment. 

3,4.2  Statistical  Analysis,  The  statistical  analysis  used  in  ANOVA  tests  begins 
with  a  sum  of  squares  approach  based  on  Equation  (3.32) 

SSt  -  SSa  +  SSb  +  SSab  +  SSe  (3.32) 

where: 

•  SSt  denotes  the  total  sum  of  squares, 

•  SSa  denotes  the  sum  of  squares  pertaining  to  variable  A, 
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•  SSb  denotes  the  sum  of  squares  pertaining  to  variable  B, 

•  SSab  denotes  the  sum  of  squares  pertaining  to  variable  interactions, 

•  SSe  denotes  the  sum  of  squares  pertaining  to  remdom  error. 

In  the  case  of  variables  A  and  B  with  n  replicates  and  the  system  response  yijk,  some 
initial  equations  must  be  defined  before  continuing. 


b  n 


yi..  —  yijk 

j=i k=i 

Vi.. 

^  Ui^ 
bn 

i  ~  1,2,, 

...,a 

a  n 

y,j,  “  X]  yijk 

i—1  k—1 

y.j. 

^  Ud:, 
an 

J-1,2, 

...,b 

n 

yij’  ~  X/ 

k^l 

Vij. 

^  IL^ 
n 

i  =  1,2,, 

...,a;  j  = 

abn 

y...  =  ^ 

y... 

^  ]L^ 
abn 

i^l  j^l  k=l 


The  dot  notation  used  above  indicates  that  the  operation  is  performed  for  all  levels  of 
the  variable  place  holder  that  the  dot  resides  in.  Formulas  for  the  sum  of  squares  in 
Equation  (3.32)  are  fisted  below. 


a  b  n 

ssT  =  T,T,Y.yii‘: 

j=l  k=l 


yl. 

abn 


^  7|2  7*2 


i=l 


bn  abn 


^  nn 


j^l 


an  abn 


For  convenience,  SSab  is  obtained  in  two  stages. 


^  V?-  7/2 

CO 

Subtotals  —  /  ^  ^  1 

“  r-f  n  abn 

1-1 


SSab  ~  SS subtotals  -  SSa  -  SSb 


SSe  -  SSt  -  SSab  -  SSa  ~  SSb 


(3.34) 

(3.35) 

(3.36) 

(3.37) 

(3.38) 

(3.39) 
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Effect 

Symbol 

Degrees  of  freedom 

A 

DOFa 

a  —  1 

B 

DOFb 

b-  1 

AB  interaction 

DOFab 

(a-l)(6-l) 

Error 

DOFe 

ab{n  —  1) 

Total 

DOFx 

abn  —  1 

Table  3.4  Degrees  of  Freedom 


Table  3.4  lists  the  degrees  of  freedom  associated  with  each  sum  of  squares.  Recalling 
that  variable  A  has  a  levels,  variable  B  has  b  levels,  and  the  experiment  has  n  replicates, 
degrees  of  freedom  indicate  the  total  number  of  variable  levels  left  that  the  variable  may 
take  on  without  remaining  in  the  same  position.  In  the  case  of  variable  interaction  AB, 
simultaneously  changing  the  levels  of  both  A  and  B  limits  the  possible  number  of  available 
variable  combinations  to  (a  —  1)(6  —  1),  thus  describing  the  number  of  degrees  of  free¬ 
dom.  In  terms  of  the  systematic  error  effect,  since  there  are  n  replicates  for  each  variable 
combination  the  system  response  has  (n  —  1)  degrees  of  freedom  for  each  variable  level 
combination.  As  a  result  for  the  entire  set  of  variable  combinations,  the  systematic  error 
has  ab{n  —  1)  degrees  of  freedom.  In  terms  of  the  total  number  of  degrees  of  freedom, 
the  two-factor  full  factorial  designed  experiment  has  abn  total  measured  system  responses 
leaving  abn  —  1  degrees  of  freedom. 

Each  sum  of  squares  divided  by  its  degrees  of  freedom  is  a  mean  square.  The  mean 
squares  applicable  to  the  ANOVA  test  pertain  to  A,  B,  AB,  and  error. 


(3.40) 

— 

(3.41) 

SSab 
=  DOF., 

(3.42) 

MS,= 

^  DOFe 

(3.43) 

Table  3.5  displays  the  resulting  ANOVA  table  and  corresponding  F  statistics.  The  F  sta¬ 
tistic  determines  whether  to  reject  or  accept  the  null  hypotheses  mentioned  in  Section  3.4.1. 
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Source  of 
Variation 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

Fo 

A  treatments 

DOFa 

MSa 

TP  _  M  US  A 

^0  “  MSe 

B  treatments 

SSb 

DOFb 

MSb 

TP.  _  MSb 

■^0  -  MSe 

Interactions 

SSab 

DOFab 

MSab 

^0  “  MSe 

Error 

SSe 

DOFe 

MSe 

Total 

SSt 

DOFt 

Table  3,5  ANOVA  Table  for  the  Two-Factor  Factorial  Experiment 


Variable  B 

Variable 

A 

15 

70 

125 

Vi.. 

1 

130  155 
74  180 

34  40 

80  75 

20  70 
82  58 

998 

2 

150  188 
159  126 

136  122 
106  115 

25  70 
58  45 

1300 

3 

138  110 
168  160 

174  120 
150  139 

96  104 
82  60 

1501 

y.j. 

1738 

1291 

770 

3799  =  y... 

Table  3.6  Experiment  Data  for  Example  ANOVA  Test 


To  determine  acceptance  or  rejection,  the  F  statistic  uses  a  confidence  interval  measure, 
a,  that  indicates  a  percentile  region  in  which  the  hypothesis  under  test  is  true.  In  the  case 
of  row  effects  due  to  A,  if  Fq  <  Fa^DOFA,DOFE^  then  accept  Hq^  and  reject  Hi,  Otherwise, 
reject  JToj  and  accept  Hi,  In  the  case  of  column  effects  due  to  jB,  if  Fq  <  Fa^nOFs^DOFE^ 
then  accept  ffo?  and  reject  Hi.  Otherwise,  reject  Hq,  and  accept  Hi.  In  the  case  of 
variable  interaction  due  to  AF,  if  Fq  <  Foc^doFab^doFe'*  accept  £fo?  and  reject  Hi. 
Otherwise,  reject  i/o?  and  accept  Hi.  Several  statistical  texts  such  as  Montgomery  and 
Box  and  Draper  include  values  for  the  F  statistic  in  tabular  form  for  confidence  intervals 
ranging  from  one  percent  to  25  percent  [5,14]. 

3.4.3  Example  Analysis  of  Variance  Test.  Using  the  ANOVA  method  just  de¬ 
scribed,  Montgomery  provides  an  example  two- factor  ANOVA  test  with  four  replicates  by 
analyzing  the  data  of  Table  3.6  [14].  The  sum  of  squares  are  computed  as  follows: 

SSt  =  (130)2  +  (155)2  +  (74)2  +  . . .  +  (60)2  _  =  77646.97  (3.44) 
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Source  of 
Variation 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

Fo 

a  =  0.05 

F  statistic 

A  treatments 

10683.72 

2 

5341.86 

7.91 

3.35 

B  treatments 

39118.72 

2 

19558.36 

28.97 

3.35 

Interactions 

9613.78 

4 

2403.44 

3.56 

2.73 

Error 

18230.75 

27 

675.21 

Total 

77646.97 

35 

Table  3.7  Example  ANOVA  Results 


(998)2  +  (1300)2  +  (1501)2  (3799)2 

~  (3)  (4)  36 

_  (1738)2  +  (1291)2  +  (770)2  (3799)2 

~  (3)(4)  36 


10683.72 

39118.72 


(3.45) 

(3.46) 


SSab 


(130  +  155  +  74  +  180)2  (34  +  40  +  80  +  75)2  (39  +  79  +  82  +  58)2 

4  +  4  +  4 

,  (150  +  188  +  159  +  126)2  (;^3g  ^22  +  106  +  115)2  ,  (25  +  70  +  58  +  45)2 

H - ;;  I - ;;  . 

4  4  4 

(138  +  110  +  168  +  160)2  (174+120+150  +  139)2  (96  +  104  +  82  +  60)2 

+  4  +  4+4 

_  (37^  _  ^9gg3  72  _  39118.72 
36 

9613.78 


SSe  =  77646.97  -  10683.72  -  39118.72  -  9613.78  =  18230.75  (3.48) 


The  analysis  of  variance  is  shown  is  Table  3.7.  Using  a  five  percent  confidence  interval,  the 
independent  variable  F  statistic  becomes 


i^o.05,2,27  =  3.35  , 


(3.49) 


and  the  variable  interaction  F  statistic  becomes 


^0.05,4,27  =  2.73  .  (3.50) 

Since  the  calculated  independent  variable  Fq  statistics  exceed  Equation  (3.49)  and  the 
calculated  variable  interaction  Fq  statistic  exceeds  Equation  (3.50),  the  effects  of  the  in- 


3-21 


dependent  variables  and  their  interactions  are  significant  within  a  95  percent  confidence 
interval.  As  a  result,  null  hypotheses  are  considered  false,  and  the  alternative  hypotheses 
are  accepted.  Each  variable  and  their  interactions  are  significant  to  the  example  experi¬ 
ment. 

3.5  Optimization  Procedures 

After  gathering,  modeling,  and  analyzing  the  system  response  data,  the  next  step 
in  this  design  methodology  implements  an  optimization  procedure  to  locate  an  optimum 
system  configuration  for  a  desired  system  response.  This  section  discusses  two  methods  of 
optimization.  The  first  method,  referred  to  here  as  the  quadratic  optimization  method,  uses 
a  technique  reserved  for  least  squares  modeling  and  is  limited  to  second  order  empirical 
models.  The  second  method  uses  the  Simplex  optimization  method  that  implements  a 
gradient  search  approach  and  applies  to  empirical  models  of  any  order. 

3.5.1  Quadratic  Optimization  Method.  Based  on  least  squares  system  modeling, 
the  quadratic  optimization  method  consists  of  two  phases:  a  linear  phase  and  a  quadratic 
phase  [14].  Using  system  models  of  the  forms: 

k 

y  =  Po  +  Yl  PiXi  +  e  (3.51) 

i—l 

and 

k  k 

y  =  Po  +  '^  /3iXi  +  ^  /3iiXi  +  ^  5^  PijXiXj  +€  (3.52) 

i—l  i=l  i  j 

> - ^ - - 

i<j 

for  the  linear  and  quadratic  phase  respectively,  the  quadratic  optimization  method  at¬ 
tempts  to  optimize  system  response  surfaces  similar  to  Figure  3.5. 

3. 5. 1.1  Sequential  Procedure.  The  quadratic  optimization  method  sequen¬ 
tially  locates  vicinities  of  system  optimum  responses  and  applies  more  elaborate  system 
models  to  continue  searching  for  the  optimum  location  [14].  Using  the  linear  equation 
to  locate  an  approximate  optimum  vicinity,  the  quadratic  optimization  method  applies  a 


3-22 


Measured  System  Response 


second-order  model  in  the  vicinity  and  continues  the  search.  Upon  locating  an  optimum 
vicinity,  the  quadratic  optimization  method  recalculates  Equation  (3.52)  over  a  smaller 
variable  region  and  proceeds  with  the  optimum  search.  Unfortunately,  for  each  new  second- 
order  model,  the  quadratic  optimization  method  requires  new  measured  system  response 
data  in  each  new  optimum  vicinity. 


3, 5 A. 2  Second-Order  Model  Analysis,  Replacing  Equation  (3,52)  with  a 
simpler  matrix  linear  algebra  equation  of  the  form 


y  f3o  +  x^b  +  x^Px  +  e 


(3.53) 


where 


x'  =  [xi,X2,xzr-- -.Xk]  , 


(3.54) 


(3.55) 
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(3.56) 


and 
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(3.57) 


[  symmetrical  f3kk  J 

allows  for  the  straightforward  location  of  optimum  vicinities.  Optimum  system  responses 
may  be  either  maximum  or  minimum  system  responses.  In  the  case  of  a  desired  maximum 
system  response,  derivatives  with  respect  to  input  variables  taken  in  the  vicinity  of  a 
maximum  system  response  quickly  indicate  its  location. 


dx 


=  b  +  2pix  =  0 


(3.58) 
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Linear  algebraic  techniques  quickly  determine  optimum  values  for  x  and  y. 


(3.59) 


^opt  —  )b 

(3.60) 

Vopt  ~  Po  i^^optb 

(3.61) 

The  accuracy  of  the  quadratic  optimization  method  depends  on  the  accuracy  of  the  second- 
order  system  model  and  the  second-order  model’s  ability  to  characterize  the  true  system 
response. 
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3.5,2  Simplex  Method,  The  simplex  optimization  method  is  a  systematic  proce¬ 
dure  of  trials  which  arrives  at  the  optimal  solution  through  the  following  iterative  steps  [16] : 

L  Determine  a  basic  feasible  solution. 

2.  Replace  the  first  basic  feasible  solution  by  another  which  improves  the  function  / 
toward  the  optimal  solution. 

3.  Continue  the  process  until  f optimal  is  obtained,  which  will  be  indicated  by  a  compu¬ 
tational  limiting  criteria,  or  until  the  computations  indicate  that  f optimal  approaches 
infinity. 

In  general,  the  function  under  consideration  contains  n  variables  and  is  described  by 

/  =  CiXl  +  C2X2  H - h  CnXn  (3.62) 

with  m  constraint  equations 

CbiiXi  +  ai2X2  + - 1-  ClinXn  <  , 

a2lXi  +  022X2  H - h  a2nXn  <  h  ^  . 

<  m<n  (3.63) 

^  ^mlXl  “h  CLm2X2  OmnXn  ^  bm  ? 

with  the  restriction  that  xj  >  0  for  all  j  [16].  In  Equation  (3.63),  xj  represents  the 
jth  design  variable,  Oij  represents  design  variable  coefficients,  and  hi  represents  the  ith 
design  limitation.  For  convenience,  the  constraint  equations  are  multiplied  by  —1  whenever 
necessary  so  that  all  bi  >  0.  Furthermore,  all  inequalities  are  changed  to  equalities  by  the 
algebraic  addition  of  new  variables  with  unit  coefficients;  that  is, 

0>ilXi  +  ai2X2  +  •  • '  +  ^in^n  i  Xji-\-i  —  bi  {%  —  1, . . .  .  (3.64) 

The  new  variables,  called  slack  variables,  do  not  contribute  to  the  value  of  the  object 
function;  that  is,  their  coefficients  Cj  are  taken  to  be  zero.  For  example,  in  the  case  of  four 
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Figure  3.6  Expansion  of  Simplex  Constraint  Equations  Towards  an  Optimal  Solution 
design  variables  with  three  constraint  equations,  n  =  4  and  m  =  3, 

n+m 

/  =  CjXj  =  Cixi  +  C2X2  +  c^xz  +  040:4  +  050:5  +  060:6  +  070:7  (3.65) 

0=1 

with  C5  =  06  =  07  =  0. 

Prom  the  m  constraint  equations,  the  system  represents  m  hyperplanes  where  typi¬ 
cally  an  improved  feasible  solution  is  a  corner  in  the  corresponding  hyperpolyhedron  [8]. 
To  locate  a  further  improved  feasible  solution,  the  constraint  equations  are  adjusted  to  shift 
the  resulting  hyperpolyhedron  in  the  direction  of  a  further  optimal  solution.  For  example, 
in  Figure  3.6  the  corner  of  the  original  hyperpolyhedron  corresponding  to  the  largest  func¬ 
tional  value  is  moved  through  the  opposite  hyperpolyhedron  face  to  a  point  with  a  lower 
functional  value.  Once  again,  a  feasible  solution  corresponding  to  an  improved  optimal 
solution  generally  resides  in  a  corner  of  the  newly  formed  hyperpolyhedron.  The  process 
continues  until  limiting  criteria  are  satisfied  and  the  final  optimal  solution  is  determined. 

3.6  Closing  Comments 

RSM  is  a  learning  and  iterative  process.  The  normal  application  of  RSM  requires 
numerous  repetitions  of  the  entire  procedure  for  accurate,  reliable  results.  Prom  the  first 
RSM  task,  constant  brainstorming  of  new  input  variables  and  consideration  of  new  experi¬ 
mental  designs  must  occur  throughout  the  process  to  guarantee  all  possible  relative  factors 
have  been  considered.  From  the  second  RSM  task,  empirical  modeling  provides  a  means 
of  predicting  system  responses  in  a  specific  region  of  interest  without  using  cumbersome 
mechanistic  models  where  the  adequacy  of  the  empirical  model  is  easily  determined  us¬ 
ing  simple  graphical  techniques.  From  the  third  RSM  task,  analysis  of  variable  effects  on 
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the  desired  system  response  using  ANOVA  tests  helps  determine  both  strong  and  weak 
variable  effects  thereby  identifying  both  necessary  system  variables  and  variables  of  lesser 
importance.  Prom  the  final  RSM  task,  optimization  of  the  system  under  study  using  an 
empirical  model  and  either  the  quadratic  optimization  method  or  the  simplex  optimiza¬ 
tion  method  provides  a  fast  means  of  computing  an  optimal  system  solution  compared 
to  tedious  mechanistic  model  testing.  The  final  comment,  RSM  is  not  intended  to  pro¬ 
vide  increased  insight  into  system  operations,  but  rather,  RSM  simply  describes  system 
performance  as  it  relates  to  various  input  variables. 
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IV,  Optimized  Design  Methodology  Development 


Using  the  properties  of  dielectric  overlays  and  the  methods  of  Response  Surface 
Methodology  with  designed  experiments  as  outlined  in  Chapters  II  and  III,  this  chapter 
develops  an  optimized  design  methodology  for  the  cavity-backed  microstrip  antenna  of 
Figure  4.1  by  implementing  a  full  factorial  designed  experiment,  cubic  spline  empirical 
modeling,  Analysis  of  Variance  (ANOVA)  tests,  and  simplex  method  optimization.  The 
design  methodology  seeks  an  optimized  antenna  design  with  minimum  average  antenna  gain 
above  the  antenna’s  bore  sight.  Section  4.1  develops  the  full  factorial  designed  experiment, 
and  Section  4.2  compares  the  adequacies  of  the  cubic  spline  empirical  model  to  those  of 
a  second  order  least  squares  empirical  model.  Section  4.3  performs  the  ANOVA  tests 
indicating  variable  effects  introduced  on  the  measured  system  response,  and  Section  4.4 
locates  an  optimized  antenna  design  through  the  use  of  the  simplex  optimization  method. 
Finally,  Section  4.5  studies  empirical  modeling  of  antenna  input  impedance  and  studies 
average  antenna  gain  empirical  modeling  close  to  a  metallic  cavity  ground  plane. 

4.1  Experiment  Development 

Prom  Chapter  III,  the  first  Response  Surface  Methodology  (RSM)  task  designs  a 
series  of  experiments  that  yields  reliable  system  response  measurements  in  a  design  region 
of  interest.  The  task  of  measuring  reliable  system  responses  requires  answers  to  the  five 
questions  of  Section  3.2.  This  section  provides  solutions  to  each  question  as  it  relates  to 
Figure  4.1  and  develops  the  needed  parameters  for  the  creation  of  a  four- variable,  five-level, 
full  factorial  designed  experiment. 

4.1.1  Input  Variables,  The  first  question  of  Section  3.2  addresses  the  issue  of 
designed  experiment  input  variables.  The  antenna  design  of  Figure  4.1  provides  several 
possibilities  of  experimental  variable  selections  upon  which  to  build  a  designed  experiment. 
A  feasible  number  of  system  response  measurements,  empirical  system  modeling,  and  time 
constraints  limit  the  number  of  designed  experiment  variables  to  four  variables  which 
include: 

1.  Antenna  operating  frequency. 
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Cavity  Dimensions: 

Length 

Width 

Height 

11.75  cm 

11.75  cm 
2.37  cm 

Patch  Dimensions: 

Length 

6.53  cm 

Width 

6.53  cm 

Patch  Location 

Centered  0.711  cm 
Above  Cavity  Bottom 

Feed  Pin  Location 

0.653  cm  Off  Center 

Substrate  Relative  Permittivity 

4.0 

Substrate  Thickness 

0.948  cm 

Table  4,1  Constant  Antenna  Design  Parameters 


2.  First  dielectric  superstrate  layer  relative  permittivity, 

3.  Second  dielectric  superstrate  layer  relative  permittivity, 

4.  First  dielectric  superstrate  layer  thickness. 

Table  4.1  lists  the  remaining  antenna  design  parameters  and  their  respective  constant 
design  values. 

Variable  Screening  Procedure.  The  second  question  of  Section  3.2  addresses 
the  importance  of  each  individual  variable  parameter  used  in  the  designed  experiment. 
Based  on  information  provided  in  Chapter  II  and  system  knowledge,  the  four  chosen  exper¬ 
iment  variables  strongly  influence  all  possible  antenna  responses  of  Figure  4.1.  Therefore, 
the  antenna  operating  frequency,  the  first  dielectric  superstrate  layer  relative  permittiv¬ 
ity,  the  second  dielectric  superstrate  layer  relative  permittivity,  and  the  first  dielectric 
superstrate  layer  thickness  are  excellent  choices  of  designed  experiment  variables. 

4^1.3  Variable  Levels.  The  third  question  of  Section  3.2  addresses  the  number  of 
levels  of  each  designed  experiment  variable.  The  number  of  levels  for  each  variable  directly 
impacts  the  ability  of  the  designed  experiment  to  capture  the  measured  system  response 
under  all  possible  variable  conditions.  Unfortunately,  a  large  number  of  variable  levels 
requires  a  greater  number  of  system  response  measurements  and  results  in  increased  time 
expenditures.  Therefore,  for  an  adequate  coverage  of  the  antenna  response  space  within  a 
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feasible  time  limit,  this  design  methodology  limits  the  variable  level  number  to  five  levels 
for  each  variable. 

4.1.3. 1  Antenna  Operating  Frequency.  The  first  designed  experiment  vari¬ 
able,  antenna  operating  frequency,  necessitated  two  initial  antenna  experiments  for  an 
adequate  selection  of  frequency  levels.  Using  a  Hybrid  Finite  Element  Method  (HFEM) 
electromagnetic  analysis  code  and  the  antenna  configuration  of  Figure  4.1,  the  first  ex¬ 
periment  set  both  dielectric  superstrate  layer  relative  permittivity  values  equal  to  1.0  and 
performed  gain  and  input  impedance  calculations  in  the  frequency  range  of  1.0  GHz  to 
1.4  GHz.  Following  the  same  procedure  as  the  first  experiment,  the  second  experiment 
set  both  dielectric  superstrate  relative  permittivities  to  8.0  and  repeated  the  gain  and  im¬ 
pedance  calculations  in  the  frequency  range  from  0.9  GHz  to  1.4  GHz.  Converting  the 
antenna  input  impedance  measurements  from  the  two  preliminary  experiments  to  antenna 
input  admittances  and  plotting  the  admittances  using  Smith  charts  normalized  to  the  ad¬ 
mittance  corresponding  to  50  ohms.  Figure  4.2  plots  admittances  from  the  first  frequency 
experiment,  and  Figure  4.3  plots  admittances  from  the  second  experiment  according  to 
their  respective  frequency  values.  Interpolating  from  both  experiments,  the  antenna  in¬ 
put  admittance  approaches  the  normalized  1.0  mho  circle  in  the  frequency  range  bounded 
by  1.2  and  1.3  GHz.  As  a  result,  using  1.3  GHz  as  an  upper  frequency  limit  led  to  the 
frequency  level  selections  consisting  of  0.9,  1.0,  1.1,  1.2,  and  1.3  GHz. 

4. 1.3.2  First  and  Second  Dielectric  Superstrate  Layer  Permittivities.  Level 
selections  of  the  second  and  third  designed  experiment  variables,  the  first  and  second  di¬ 
electric  superstrate  layer  relative  permittivities  from  Figure  4.1,  came  about  from  practical 
antenna  design  experience.  Common  dielectric  materials  with  relative  permittivity  values 
in  the  range  of  1.0  to  8,0  are  readily  available  or  at  least  manufacturable.  Therefore,  the 
accessibility  of  common  dielectric  materials  led  to  the  selections  of  1.0,  2.0,  4.0,  6.0,  and 
8.0  for  dielectric  relative  permittivity  values  used  in  the  designed  experiment. 

4. 1.3.3  First  Dielectric  Superstrate  Layer  Thickness.  Using  an  HFEM  an¬ 
tenna  analysis  code  and  the  spatial  partitioning  property  of  HFEM  codes  within  cavity 


4-4 


Figure  4.2 


Figure  4.3 


Smith  Chart  Analysis  of  Input  Admittance  Values  With  First  and  Second 
Superstrate  Dielectric  Relative  Permittivities  Equal  to  1.0  {eri  =  er2  =  1-0) 


Smith  Chart  Analysis  of  Input  Admittance  Values  With  First  and  Second 
Superstrate  Dielectric  Relative  Permittivities  Equal  to  8.0  =  er2  —  8.0) 


4-5 


Frequency  (GHz) 

0.9,  1.0,  1.1 
1.2,  and  1.3 

First  Layer  Relative 
Permittivity 

1,  2,  4, 
6,  and  8 

Second  Layer  Relative 
Permittivity 

1,  2,  4, 
6,  and  8 

First  Layer 

Thickness  (cm) 

0.237,  0.474, 
0.711,  0.948, 
and  1.185 

Table  4.2  Variable  Antenna  Design  Parameters 

dimensions,  the  level  selections  for  the  fourth  designed  experiment  variable,  the  first  dielec¬ 
tric  superstrate  layer  thickness,  was  limited  to  HFEM  nodal  dimensions.  The  first  dielectric 
superstrate  layer  thickness  levels  include  0.237,  0.474,  0.711,  0.948,  and  1.185  centimeters. 
Worth  noting,  due  to  constant  cavity  dimensions  and  a  constant  substrate  depth,  variations 
in  the  first  superstrate  thickness  directly  impact  the  thickness  of  the  second  superstrate 
layer  causing  the  second  superstrate  layer  to  compensate  for  the  remaining  cavity  space. 
Table  4.2  lists  the  designed  experiment  variables  and  their  corresponding  variable  levels. 

System  Response,  In  this  development  of  an  optimized  design  methodology 
for  the  antenna  in  Figure  4.1,  the  desired  system  response  is  the  average  antenna  gain  over 
the  angular  region  described  in  Equations  (4.1)  and  (4.2). 


0^  <  6^  <  20° 


(4.1) 


0°  <  (^  <  360°  (4.2) 

Using  gain  calculations  in  five  degree  increments.  Equation  (4.3)  defines  average  antenna 
gain. 


/  /  gairif^o  \  /  gain^o  \ 

/  gain^f^o  \  \ 

averagegain  =  20  log 

10 V  /  +  lOV  ^  +  .. 

..  +  10V  / 

5 

J 

4 A. 5  Experimental  Design.  The  fourth  question  of  Section  3.2  addresses  the  type 
of  experimental  design  used  to  gather  measured  system  response  data.  The  electromag¬ 
netic  properties  of  the  antenna  design  in  Figure  4.1  and  antenna  radiation  pattern  lobing 
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Test 

Layer  1 
(cm) 

Crl 

€r2 

Preq 

(GHz) 

Gain 

(dB) 

1 

0.237 

1 

1 

0.9 

yi 

2 

0.237 

1 

1 

1.0 

2/2 

3 

0.237 

1 

1 

1.1 

ys 

4 

0.237 

1 

1 

1.2 

y4 

5 

0.237 

1 

1 

1.3 

y5 

6 

0.237 

1 

2 

0.9 

ye 

625 

1.185 

8 

8 

1.3 

2/625 

Table  4.3  Fall  Factorial  Design 


necessitate  a  full  factorial  designed  experiment.  The  full  factorial  experiment  meetsures 
the  average  antenna  gain  response  at  each  variable  level  combination.  Since  five  variable 
levels  for  each  of  the  four  variables  were  selected,  the  full  factorial  experiment  determines 
the  average  antenna  gain  at  625  (5^)  total  experimental  locations.  Table  4.3  lists  the  full 
factorial  experiment  used  in  this  experiment. 

4^1*6  Hybrid  Finite  Element  Method  Measurements.  The  fifth  question  of  Sec¬ 
tion  3.2  addresses  the  means  of  measuring  the  system  response.  In  this  optimized  design 
methodology  development,  an  HFEM  electromagnetic  analysis  code  directly  suited  for 
rectangular  cavity  aperture  antennas  is  used  that  uses  brick  segments  to  facilitate  field 
solutions  in  the  cavity  region  and  integral  equations  to  solve  for  field  quantities  in  the 
unbounded  external  region. 

4.2  Average  Antenna  Gain  Empirical  Modeling 

From  Chapter  III,  the  second  RSM  task  determines  an  empirical  model  that  accu¬ 
rately  describes  the  measured  system  response,  the  average  antenna  gain.  Several  possible 
empirical  models  exist  of  which  cubic  splines  and  least  squares  represent  two  common 
techniques.  Section  4.2.1  addresses  the  application  of  cubic  spline  empirical  modeling 
to  this  multi-dimensional  space  problem,  and  Section  4.2.2  provides  a  second  order  least 
squares  empirical  model  applicable  to  the  optimization  technique  discussed  in  Section  3.5.1. 
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Section  4.2.3  compares  the  effectiveness  of  the  cubic  spline  empirical  model  with  the  effec¬ 
tiveness  of  the  least  squares  empirical  model. 

4.2.1  Cubic  Splines  Empirical  Modeling.  The  application  of  cubic  splines  in  a 
multi-dimensional  space  requires  an  extension  of  the  cubic  spline  derivation  of  Section  3.3.3. 
Utilizing  a  small  portion  of  the  Wilson  [20]  cubic  spline  algorithm  that  extends  the  theory 
of  Section  3.3.3  to  multi-dimensional  space,  the  cubic  spline  empirical  modeling  algorithm 
performs  the  following  sequential  steps: 

1.  Use  the  initial  625  (5^)  test  configurations  of  the  full  factorial  designed  experiment 
and  their  corresponding  gain  values  and  creates  cubic  splines  in  the  direction  of  the 
first  input  variable,  the  antenna  operating  frequency,  at  all  possible  level  combina¬ 
tions  of  the  remaining  three  input  variables. 

2.  Set  the  antenna  operating  frequency  variable  to  its  desired  input  value. 

3.  Determine  average  antenna  gain  values  in  terms  of  the  remaining  three  input  variables 
and  their  corresponding  five  levels  at  the  desired  antenna  operating  frequency  for  a 
total  of  125  (5^)  new  average  antenna  gain  values. 

4.  Use  the  resulting  125  possible  test  configurations  and  their  corresponding  gain  values 
and  creates  cubic  splines  in  the  direction  of  the  second  input  variable,  the  first  dielec¬ 
tric  superstrate  layer  permittivity,  at  all  possible  level  combinations  of  the  remaining 
two  input  variables. 

5.  Set  the  first  dielectric  superstrate  layer  permittivity  value  to  its  desired  input  value. 

6.  Determine  average  antenna  gain  values  in  terms  of  the  remaining  two  input  variables 
and  their  corresponding  five  levels  at  the  desired  input  values  of  both  the  antenna 
operating  frequency  and  the  first  dielectric  superstrate  layer  permittivity  for  a  total 
of  25  (5^)  calculations. 

7.  Use  the  resulting  25  possible  test  configurations  and  their  corresponding  gain  values 
and  creates  cubic  splines  in  the  direction  of  the  third  input  variable,  the  second 
dielectric  superstrate  layer  permittivity,  at  all  possible  levels  of  the  remaining  input 
variable. 
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8.  Set  the  second  dielectric  superstrata  layer  permittivity  value  to  its  desired  input 
value. 

9.  Determine  average  antenna  gain  values  in  terms  of  the  remaining  input  variable,  the 
first  dielectric  superstrata  layer  thickness,  and  its  corresponding  five  levels  at  the 
desired  input  values  of  antenna  operating  frequency,  first  dielectric  superstrata  layer 
permittivity,  and  second  dielectric  superstrata  layer  permittivity  for  a  total  of  five 
calculations. 

10.  Use  the  resulting  five  possible  test  configurations  and  their  corresponding  gain  values 
and  creates  cubic  splines  in  the  direction  of  the  final  input  variable,  first  dielectric 
superstrata  layer  thickness. 

11.  Set  the  first  dielectric  superstrata  layer  thickness  to  its  desired  input  level  and  cal¬ 
culates  the  final  average  antenna  gain  value  that  now  represents  the  desired  levels  of 
all  four  input  variables. 

The  multi-dimensional  cubic  spline  empirical  model  provides  a  computationally  efficient 
model  which  can  be  optimized  using  the  simplex  optimization  method. 

4^2,2  Least  Squares  Empirical  Modeling,  Using  the  average  antenna  gain  mea¬ 
surements  obtained  from  the  full  factorial  designed  experiment  of  Section  4.1.5  and  a  least 
squares  software  routine  from  Matlab,,  Equation  (4.4)  defines  a  second-order  quadratic 
equation  describing  the  antenna  of  Figure  4.1  for  the  variable  ranges  of  Table  4.2. 

Vgain  =  -11.2938(/r'eg)2  +  0.0018(e^2)^  (4.4) 

+0.0004(eri)^  -  0.0050{thickif 
-0A264{freq){er2)  -  0.S090{freq){€ri) 
+0.0479{freq){thicki)  —  0.0067(€r2)(€ri) 

+0.0105(er2)(</i*cfei)  —  0.0105(eri)(i/i.icfei) 

+27.6203(/req-)  +  0.4040(6^2)  +  0.3682(eri) 

-0.0189{thicki)  -  14.6896 
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Frequency  (GHz) 

0.95,  1.15,  1.25 

First  Layer 
Permittivity 

1.5,  5,  7 

Second  Layer 
Permittivity 

1.5,  5,  7, 

First  Layer 
Thickness  (cm) 

0.237,  0.711,  1.185 

Table  4.4  Model  Adequacy  Test  Variable  Levels 


Test 

Layer  1 
(cm) 

tr\ 

€r2 

Preq 

(GHz) 

Gain 

(dB) 

1 

0.237 

1.5 

1.5 

0.95 

yi 

2 

0.237 

1.5 

1.5 

1.15 

2/2 

3 

0.237 

1.5 

1.5 

1.25 

ys 

4 

0.237 

1.5 

5 

0.95 

Vi 

81 

1.185 

7 

7 

1.25 

yu 

Table  4.5  Model  Adequacy  Test  Cases 


^.2.3  Average  Antenna  Gain  Model  Adequacy,  Using  a  set  of  81  test  cases  from 
Table  4.4  and  Table  4.5,  the  techniques  of  Section  3.3.4  pictorially  describe  the  empirical 
modeling  adequacies  of  the  cubic  spline  and  least  squares  models. 

4^2,3. 1  Normal  Probability  Plot  of  Residuals.  Using  the  normal  probability 
plot  technique  of  Section  3.3.4. 1,  Figure  4.4  indicates  that  the  cubic  spline  empirical  model 
accuracy  far  exceeds  that  of  the  least  squares  empirical  model.  In  Figure  4.4,  since  the 
cumulative  probabilities  of  the  cubic  spline  model  residuals  approximate  a  straight  line 
better  than  the  cumulative  probabilities  of  the  least  squares  model  residuals,  the  normally 
distributed  systematic  error  in  the  cubic  spline  model  is  less  than  the  systematic  error  in 
the  least  squares  model.  Furthermore,  from  the  cubic  spline  empirical  model,  the  average 
antenna  gain  residuals  vary  over  a  net  range  of  1.2  dB,  while  residuals  from  the  least 
squares  empirical  model  vary  over  a  net  range  of  approximately  4.5  dB,  a  significantly 
larger  gain  region  than  the  cubic  spline  model.  Table  4.6  lists  the  residual  mean  and  the 
residual  standard  deviation  for  both  empirical  models. 
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Least  Squares  Probability  Cubic  Spline  Probability 


Normal  Probability  Plot 


Figure  4.4  Normal  Probability  Plots  of  Cubic  Spline  and  Least  Squares  System  Response 
Models 


Cubic  Spline: 

Residual  Mean 

Residual  Standard  Deviation 

0.0357  dB 
0.1592  dB 

Least  Squares: 

Residual  Mean 

Residual  Standard  Deviation 

-0.0579  dB 
0.8240  dB 

Table  4.6  Mean  and  Standard  Deviation  Values  of  Cubic  Spline  and  Least  Squares  Resid¬ 
uals 
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Cubic  Spline  Fitted  Values  (dB) 


Figure  4.5  Cubic  Spline  and  Least  Squares  Residuals  Versus  Fitted  Values 

4^2. 3,2  Residuals  Versus  Fitted  Value  Plots,  Using  the  graphical  tech- 
niques  of  Section  3.3.4.2,  the  plots  of  empirical  model  residuals  versus  the  empirical  model 
predicted  average  antenna  gain  values  in  Figxire  4.5  indicate  that  the  cubic  spline  model 
residuals  vary  much  less  than  those  from  the  least  squares  model.  As  indicated  in  the  nor¬ 
mal  probability  plots  in  Figure  4.4,  the  residual  versus  fitted  value  plots  indicate  that  the 
cubic  spline  model  residuals  vary  over  a  net  range  of  1.2  dB,  and  the  least  squares  model 
residuals  vary  over  a  net  range  of  4.5  dB.  Since  both  plots  in  Figure  4.5  are  structure-less, 
both  empirical  models  indicate  a  normal  random  distribution  of  model  residuals. 
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Antenna  Frequency  (GHz) 


Figure  4.6  Cubic  Spline  and  Least  Squares  Residuals  Versus  Frequency  Values 


^,2.3.3  Residuals  Versus  Variable  Parameter  Plots.  Using  the  graphical 
technique  of  Section  3. 3.4.3,  the  plots  of  residuals  versus  input  variable  parameters  in 
Figures  4.6,  4.7,  4.8,  and  4.9  again  indicate  that  the  cubic  spline  empirical  model  out^ 
performs  the  least  squares  empirical  model.  As  before,  residuals  from  the  cubic  spline 
empirical  model  vary  over  an  approximate  net  range  of  1.2  dB,  and  residuals  from  the 
least  squares  empirical  model  vary  over  an  approximate  net  range  of  4.5  dB. 

From  the  model  adequacy  plots,  cubic  spline  empirical  modeling  far  exceeds  the  least 
squares  performance  and  will  be  implemented  throughout  the  remainder  of  this  experiment. 
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Residuals  (dB) 


J^,3  Analysis  of  Variance 

Prom  Chapter  III,  the  third  RSM  task  analyzes  the  measured  system  response  in 
relation  to  the  effects  introduced  by  the  designed  experiment  input  variables  and  their 
interactions.  Section  4.3.1  implements  the  ANOVA  techniques  of  Section  3.4  to  analyze 
the  effects  introduced  by  the  antenna  operating  frequency  on  the  average  antenna  gain 
by  using  average  antenna  gain  predictions  from  the  cubic  spline  empirical  model.  After 
completing  the  cubic  spline  based  ANOVA  test  for  frequency  effects,  Section  4.3.1  then 
validates  the  frequency  test  by  reaccomplishing  the  frequency  ANOVA  test  using  average 
antenna  gain  values  from  the  HFEM  code  and  compares  findings.  The  purpose  of  repeating 
the  frequency  ANOVA  test  using  HFEM  results  is  to  demonstrate  the  reliability  of  the  cubic 
spline  empirical  model  for  ANOVA  testing  thereby  minimizing  the  computational  effort 
and  time  expenditures  brought  on  by  using  HFEM  results.  After  proving  the  effectiveness 
of  cubic  spline  empirical  modeling  for  ANOVA  tests.  Section  4.3.2  accomplishes  an  ANOVA 
test  studying  the  effects  of  all  four  input  variables  as  well  as  their  interactions  using  the 
cubic  spline  empirical  model. 

4-3. 1  Frequency  Effects  ANOVA  Test  Applying  the  technique  of  Section  3.4,  this 
section  demonstrates  the  usefulness  of  empirical  modeling  on  ANOVA  testing  by  comparing 
results  from  a  cubic  spline  empirical  model  based  frequency  ANOVA  test  to  results  from 
a  HFEM  based  frequency  ANOVA  test.  Using  the  test  hypotheses  of  Section  3.4.1,  the 
frequency  effects  ANOVA  test  attempts  to  prove  one  of  two  possible  hypotheses  true. 
The  null  hypothesis,  i?o?  specifies  that  none  of  the  frequency  levels  used  in  the  ANOVA 
test  effect  the  average  antenna  gain,  while  the  alternate  hypothesis,  Hi^  specifies  that  at 
least  one  of  the  frequency  levels  effect  the  antenna  gain.  Using  identical  normal  randomly 
distributed  sets  of  ten  replicates  for  each  frequency  level  for  the  variable  values  given  in 
Table  4.7,  Tables  4.8  and  4.9  summarize  the  results  of  the  two  frequency  ANOVA  tests.  In 
each  case,  the  calculated  Fq  statistic  exceeds  the  required  F  statistic  with  the  one  percent 
confidence  level  and  proves  the  null  hypotheses  false.  These  frequency  ANOVA  tests  results 
demonstrate  that  the  cubic  spline  empirical  model  adequately  predicts  ANOVA  results  and 
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Frequency  (GHz) 

0.95,  1.05, 
1.15,  and  1.25 

First  Layer  Relative  Permittivity 

4 

Second  Layer  Relative  Permittivity 

4 

First  Layer  Thickness  (cm) 

0.711 

Replicates  Per  Frequency  Level 

10 

Table  4.7  Frequency  ANOVA  Test  Variable  Levels 


Factor 

Sum  of 
Squsires 

Degrees  of 
Freedom 

Mean 

Square 

Fo 

a  =  0.01 

F  Statistic 

Frequency 

23.79 

3 

7.93 

14.96 

4.31 

Error 

18.94 

36 

0.53 

Total 

42.73 

39 

Table  4.8  Cubic  Spline  Empirical  Model  Frequency  ANOVA  Test  Results 

is  well  suited  to  perform  a  full  ANOVA  test  consisting  of  the  four  variables  used  in  the  full 

factorial  designed  experiment  and  their  possible  interactions. 

\ 

4-3,2  Full  Effects  ANOVA  Test  Since  Section  4.3.1  proved  the  usefulness  of 
the  cubic  spline  empirical  model  for  ANOVA  testing,  this  section  uses  the  cubic  spline 
gain  predictions  to  extend  the  ANOVA  test  in  Section  4.3.1  to  include  all  foxur  variable 
parameters  used  in  the  full  factorial  designed  experiment.  The  full  ANOVA  test  includes 
interactions  among  the  variables  and  studies  each  factor’s  effect  on  the  average  antenna 
gain.  Using  100  normal  randomly  distributed  replicates  of  the  variable  combinations  listed 
in  Table  4.10,  the  ANOVA  test  proves  that  each  variable  used  in  the  full  factorial  design 
strongly  influences  the  average  antenna  gain  for  the  antenna  design  of  Figure  4.1  in  the 
angular  region  described  by  Equations  (4.1)  and  (4.2). 


Factor 

Sum  of 
Squares 

Degrees  of 
Freedom 

Mean 

Square 

Fo 

a  =  0.01 

F  Statistic 

Frequency 

24.89 

3 

8.30 

14.07 

4.31 

Error 

21.19 

36 

0.59 

Total 

46.08 

39 

Table  4.9  HFEM  Frequency  ANOVA  Test  Results 


4-18 


Frequency  (GHz) 

0.95,  1.05, 
1.15,  and  1.25 

First  Layer  Relative 
Permittivity 

1.5,  3, 
5,  and  7 

Second  Layer  Relative 
Permittivity 

1.5,  3, 
5,  and  7 

First  Layer 

Thickness  (cm) 

0.356,  0.593, 
0.830,  and  1.067 

Replicates  Per  Combination 

100 

Table  4.10  Full  ANOVA  Test  Variable  Levels 


Full  ANOVA  Test  Hypotheses.  Extending  the  test  hypotheses  of 
Section  3.4.1  to  four  variables,  the  null  hypothesis,  Hq^  for  each  variable  parameter  pre¬ 
dicts  that  none  of  the  variable  levels  effect  the  average  antenna  gain,  while  the  alternate 
hypothesis,  Ifi,  of  each  variable  predicts  that  at  least  one  variable  level  effects  the  average 
antenna  gain.  In  terms  of  variable  interaction,  for  each  variable  interaction  combination 
the  null  hypothesis  predicts  that  none  of  the  variable  levels  included  in  the  interaction 
combination  effect  the  average  antenna  gain,  while  the  alternate  hypothesis  predicts  that 
at  least  one  variable  level  included  in  the  interaction  combination  effects  the  measured 
average  antenna  gain.  In  other  words,  the  ANOVA  test  determines  the  sensitivity  of  the 
antenna  gain  to  changes  in  the  antenna  operating  frequency,  the  first  and  second  dielectric 
superstrate  relative  permittivities,  and  the  first  dielectric  superstrate  thickness  as  well  as 
their  interactions. 

4*3. 2. 2  Full  ANOVA  Statistical  Analysis.  Applying  the  statistical  tech¬ 
niques  of  Section  3.4.2  to  the  four  experimental  variables.  Table  4.11  summarizes  the 
results  of  the  full  ANOVA  test.  In  each  case,  the  calculated  Fq  statistic  exceeds  the  re¬ 
quired  F  statistic  with  the  confidence  level  of  1  percent  and  proves  the  null  hypotheses 
false.  Each  variable  factor  and  variable  interaction  factor  significantly  affects  the  average 
antenna  gain.  Prom  Table  4.11,  frequency  has  the  greatest  influence  on  the  average  an¬ 
tenna  gain  for  the  cavity-backed  antenna  of  Figure  4.1.  Since  cavities  resonate  at  specific 
frequencies  according  to  their  dimensions  in  frequency  wavelengths,  the  frequency  effects 
for  the  full  ANOVA  test  follow  engineering  expectations.  The  relative  permittivity  of  the 
second  superstrate  layer  has  the  second  largest  influence  on  the  average  antenna  gain  of 
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Factor 

Number 

Pactor 

Sum  of 
Squares 

Degrees  of 
Preedom 

Mean 

Square 

Fo 

a  =  0.01 

F  Statistic 

1 

freq 

8018.94 

3 

2672.98 

4135.62 

3.8 

2 

Crl 

1892.91 

3 

630.97 

976.23 

3.8 

3 

€r2 

4317.31 

3 

1439.10 

2226.58 

3.8 

4 

thicki 

50.78 

3 

16.93 

26.19 

3.8 

5 

freqeri 

5147.94 

9 

571.99 

884.99 

2.4 

6 

freqer2 

10821.47 

9 

1202.39 

1860.33 

2.4 

7 

freqthicki 

121.56 

9 

13.51 

20.90 

2.4 

8 

^rl^r2 

1277.47 

9 

141.94 

219.61 

2.4 

9 

Crithicki 

320.75 

9 

35.64 

55.14 

2.4 

10 

er^thicki 

566.41 

9 

62.93 

97.37 

2.4 

11 

freqeri€r2 

1286.12 

27 

47.63 

73.70 

1.75 

12 

freqerithicki 

871.59 

27 

32.28 

49.95 

1.75 

13 

freqer2thicki 

900.25 

27 

33.34 

51.59 

1.75 

14 

ej*'\_e’p2th/tck'ji^ 

153.44 

27 

5.68 

8.79 

1.75 

15 

freqcr  i  er2thicki 

308.41 

81 

3.81 

5.89 

1.4 

Error 

14981.94 

23180 

0.65 

Total 

51037.28 

23435 

Table  4.11  Pull  ANOVA  Test  Results 


this  particular  antenna  configuration.  It  stands  to  reason  due  to  its  close  proximity  to  the 
antenna  radiating  aperture  that  the  second  dielectric  superstate  layer’s  relative  permit¬ 
tivity  has  a  stronger  influence  than  the  first  dielectric  superstrate.  In  fact,  the  variable 
interaction  created  by  the  interaction  between  anteima  operating  frequency  and  the  second 
dielectric  superstrate  relative  permittivity  has  the  third  largest  influence  on  the  average 
antenna  gain.  In  terms  of  manufacturing  tolerances  and  variable  sensitivities,  this  full 
ANOVA  test  indicates  that  tolerances  relating  to  frequency  and  the  relative  permittivity 
of  the  second  dielectric  superstrate  are  the  most  critical  and  require  tight  manufacturing 
control.  The  remaining  factors  studied  in  this  full  ANOVA  test  are  easily  ranked  accord¬ 
ing  to  their  impact  on  the  average  antenna  gain  where  each  factor  strongly  influences  the 
average  antenna  gain. 


Optimization  Procedures 

Prom  Chapter  III,  the  fourth  RSM  task  searches  for  the  optimal  levels  of  the  input 
variables  that  produce  the  desired  optimized  system  response.  In  this  experiment,  the 


4-20 


Frequency 

1.3  GHz 

First  Layer  Relative  Permittivity 

6.0 

Second  Layer  Relative  Permittivity 

8.0 

First  Layer  Thickness 

0.474  cm 

Average  Gain  (HFEM) 

-5.779  dB 

Table  4.12  Initial  Optimized  Location  Guess 


Frequency 

1.3  GHz 

First  Layer  Relative  Permittivity 

6.696 

Second  Layer  Realtive  Permittivity 

8.0 

First  Layer  Thickness 

0.512  cm 

Elapsed  Optimization  Time 

1  min 

Simplex  Iterations 

180 

Average  Gain  (Empirical  Model) 

-5.890  dB 

Average  Gain  (HFEM) 

-6.206  dB 

Input  Impedance  (HFEM) 

48.080  -h  i42.937 

Gain  Accuracy 

3.8  percent 

Gain  Improvement  Over  Initial  Guess 

4.8  percent 

HFEM  Optimization  Validation  Period 

2  hrs,  59  mins 

Table  4.13  Optimized  Antenna  Design 


desired  response  is  a  minimum  average  antenna  gain  over  the  angular  region  described 
by  Equations  (4.1)  and  (4,2).  Using  the  simplex  optimization  method  of  Section  3.5.2, 
the  cubic  spline  average  antenna  gain  empirical  model,  the  initial  optimized  design  guess 
given  in  Table  4.12,  and  the  optimization  flow  chart  of  Figure  4.10,  the  resulting  optimized 
antenna  design  is  listed  in  Table  4.13.  In  terms  of  field  quantities,  the  average  antenna 
gain  of  the  optimized  design  compares  within  3.8  percent  of  the  validated  HFEM  optimized 
design  configuration  solution  with  a  net  average  antenna  gain  reduction  of  4.8  percent  over 
the  initial  guess  design.  The  resulting  antenna  radiation  pattern  of  Figure  4.11  indicates  a 
null  located  above  the  antenna  bore  sight  in  the  desired  angular  region  thus  demonstrating 
that  this  is  a  feasible  optimized  design  methodology. 

Figures  4.12,  4.13,  4.14,  4.15,  4.16,  and  4.17  represent  the  average  antenna  gain 
response  of  the  optimized  antenna  design  as  variable  pairs  approach  their  optimized  values. 
For  example,  in  Figure  4.12  as  the  first  layer  relative  permittivity  approaches  6,696,  the 
permittivity  optimized  value,  and  as  the  antenna  operating  frequency  approaches  1.3  GHz, 
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Figure  4.10  Optimization  Process  Flow  Chart 
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Linear  "Y"  Directed  Gain 


Linear  "X"  Directed  Gain 


Figure  4.11  Optimized  Antenna  Pattern 


First  layer  permittivity 

Figure  4.12  First  Layer  Relative  Permittivity  and  Frequency  Effects  Near  the  Optimized 
Location 


the  frequency  optimized  value,  the  average  antenna  gain  approaches  a  minimum  value  near 
-6.0  dB. 

Besides  the  location  of  an  optimized  antenna  design,  Table  4.13  yields  some  inter¬ 
esting  points.  Using  the  validated  cubic  spline  empirical  model  to  locate  an  optimized 
antenna  design  for  the  antenna  configuration  of  Figure  4.1,  each  simplex  method  iteration 
took  approximately  0.33  seconds  to  complete.  When  validating  the  optimized  design  using 
an  HFEM  code,  the  single  antenna  design  configuration  took  approximately  three  hours 
to  calculate.  Assuming  each  simplex  iteration  would  take  as  long,  it  would  have  taken 
23  days  for  the  HFEM  code  to  come  up  with  the  optimized  design  that  the  RSM  method 
found  in  only  one  minute.  This  is  the  power  of  RSM  and  designed  experiments.  The  design 
methodology  used  in  this  experiment  is  accurate  and  a  tremendous  time  saving  tool. 


4^.5  Additional  System  Response  Modeling  Possibilities 

In  addition  to  the  average  antenna  gain  system  response  in  the  angular  region  de¬ 
scribed  by  Equations  (4.1)  and  (4.2),  other  possible  desired  system  responses  exist  for  the 
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Second  layer  permittivity  Frequency  (GHz) 


Figure  4.13  Second  Layer  Relative  Permittivity  and  Frequency  Effects  Near  the  Opti¬ 
mized  Location 


1 .32  2  First  layer  thickness  (m) 

Frequency  (GHz) 


Figure  4.14  First  Layer  Thickness  and  Frequency  Effects  Near  the  Optimized  Location 
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Figure  4.15  First  and  Second  Layer  Relative  Permittivity  Effects  Near  the  Optimized 
Location 


Second  layer  permittivity  First  layer  thickness  (m) 


Figure  4.16  First  Layer  Thickness  and  Second  Layer  Relative  Permittivity  Effects  Near 
the  Optimized  Location 
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Figure  4.17  First  Layer  Relative  Permittivity  and  Thickness  Effects  Near  the  Optimized 
Location 


antenna  design  of  Figure  4.1  such  as  antenna  input  impedance  or  average  antenna  gain 
near  the  cavity  ground  plane.  Using  the  full  factorial  designed  experiment  of  Section  4.1.5 
and  the  model  adequacy  tests  in  Section  4.2.3,  Section  4,5.1  analyzes  the  ability  of  the  cu¬ 
bic  spline  empirical  modeling  technique  to  adequately  describe  the  performance  of  antenna 
input  impedance  within  the  input  variable  ranges  of  Table  4.2.  Likewise,  Section  4.5.2  an¬ 
alyzes  the  ability  of  the  cubic  spline  empirical  modeling  technique  to  adequately  describe 
the  average  antenna  gain  in  the  angular  region  described  by  Equations  (4.5)  and  (4.6) 
within  the  input  variable  ranges  of  Table  4.2. 


70°  <0<  85° 

(4.5) 

0°  <4><  360° 

(4.6) 

4.5.1  Input  Impedance  Modeling.  At  first  glance,  the  impedance  modeling  ade¬ 
quacy  plots  of  Figures  4.18,  4.19,  4.20,  4.21,  4.22,  and  4.23  indicate  that  the  cubic  spline 
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Real  Impedance  Residual  Mean 

Real  Impedance  Residual  Standard  Deviation 

-0.0941  ohms 
4.4942  ohms 

Imaginary  Impedance  Residual  Mean 

Imaginary  Impedance  Residual  Standard  Deviation 

0.9459  ohms 
6.6533  ohms 

Table  4.14  Real  and  Imaginary  Impedance  Residual  Means  and  Standard  Deviations 

empirical  model  poorly  characterizes  the  antenna’s  input  impedance  performance.  In  fact, 
from  the  normal  probability  plots  both  the  real  and  imaginary  impedance  residuals  vary 
over  a  large  range  of  impedances.  Likewise,  from  the  residuals  versus  fitted  value  plots 
and  the  residuals  versus  individual  variable  plots  both  the  real  and  imaginary  impedance 
residuals  vary  over  an  excessively  large  impedance  range.  Furthermore,  Table  4.14  lists  the 
mean  impedance  residuals  along  with  corresponding  residual  standard  deviations  for  both 
the  real  and  imaginary  impedance  components  indicating  the  poor  modeling  performance 
of  the  cubic  spline  empirical  model.  Despite  the  high  modeling  variability,  one  encour¬ 
aging  prospect  remains.  The  impedance  empirical  models  predict  impedance  trends,  and 
if  employed  correctly  in  a  weighted  or  combined  optimization  procedure,  the  impedance 
empirical  models  may  provide  a  useful  optimization  tool.  Chapter  V  employs  a  combined 
optimization  procedure  composed  of  both  average  antenna  gain  and  antenna  impedance 
empirical  models. 


4-5.2  Average  Antenna  Gain  Modeling  Near  the  Cavity  Ground  Plane.  In  a 
similar  situation,  the  modeling  adequacy  plots  of  Figures  4.24,  4.25,  4.26,  4.27,  4.28,  and 
4.29  give  little  indication  that  the  four  variable  cubic  spline  empirical  model  accurately 
describes  the  average  antenna  gain  near  the  cavity  ground  plane.  Using  Equations  (4.5) 
and  (4.6),  the  average  antenna  gain  near  the  antenna  ground  plane  is  described  by 


averagegain  =  20  log 


f  lOV  ^°  /  +  lOV  / 


+  10 

4 


I  20  ) 


(4.7) 


The  antenna  gain  normal  probability  plot  of  Figure  4.24  shows  that  the  average  antenna 
gain  residuals  vary  over  a  total  range  of  approximately  8.2  dB  which  is  repeated  in  the  resid¬ 
uals  versus  fitted  values  plot  and  the  residuals  versus  individual  variable  plots.  Table  4.15 
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Probability  Probability 


Normal  Probability  Plot 


Normal  Probability  Plot 


Figure  4.18 


Normal  Probability  Plots  of  Cubic  Spline  Input  Impedance  Residuals 
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Imaginary  Residuals  (ohms)  Real  Residuals  (ohms) 


Antenna  Frequency  (GHz) 


Figure  4.20  Cubic  Spline  Impedance  Residuals  Versus  Frequency  Values 
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Second  Layer  Relative  Permittivity 

Cubic  Spline  Impedance  Residuals  Versus  Second  Layer  Relative  Permittiv 
ity  Values 


Imaginary  Residuals  (ohms)  Real  Residuals  (ohms) 
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Cubic  Spline  Model 
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Figure  4.24  Cubic  Spline  Normal  Probability  Plot  of  Residuals  Near  Ground  Plane 

Residual  Mean  0.0357  dB 

Residual  Standard  Deviation  0.1592  dB 

Table  4.15  Average  Antenna  Gain  Cubic  Spline  Model  Residual  Mean  And  Residual 
Standard  Deviation  Values  Near  Ground  Plane 

lists  the  residual  mean  and  residual  standard  deviation  values  of  the  average  antenna  gain 
empirical  model  near  the  cavity  ground  plane. 

Empirical  modeling  near  the  antenna  ground  plane  requires  smaller  variable  steps  in 
the  full  factorial  designed  experiment  to  adequately  predict  average  antenna  gain.  Lobing 
effects  due  to  antenna  pattern  side  lobes  and  metallic  cavity  interaction  make  empirical 
modeling  more  strenuous.  In  addition,  using  a  narrower  angular  region  such  as 

85°  <  0  <  90°  (4.8) 

0°  <  (/>  <  360°  (4.9) 

with  average  antenna  gain  measurements  performed  at  one  degree  increments  could  possi¬ 
bly  provide  a  better  opportunity  to  accurately  characterize  the  antenna  performance  near 
the  ground  plane. 
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Cubic  Spline  Residuals  (dB) 


Cubic  Spline  Fitted  Values  (dB) 


Figure  4.25 


Cubic  Splines  Residuals  Near  Ground  Plane  Versus  Fitted  Values 


Figure  4.26  Cubic  Splines  Residuals  Near  Ground  Plane  Versus  Frequency  Values 
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First  Layer  Relative  Permittivity 


Figure  4.27  Cubic  Spline  Residuals  Near  Ground  Plane  Versus  First  Layer  Relative  Per¬ 
mittivity  Values 
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Figure  4,28  Cubic  Spline  Residuals  Near  Ground  Plane  Versus  Second  Layer  Relative 
Permittivity  Values 
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First  Dielectric  Layer  Thickness  (m)  ^  ^  q-3 

Figure  4.29  Cubic  Spline  Residuals  Near  Ground  Plane  Versus  First  Layer  Thickness 
Values 

4-6  Closing  Comments 

Developed  here,  an  RSM  based  optimized  antenna  design  methodology  composed 
of  a  full  factorial  designed  experiment,  a  cubic  spline  empirical  model,  empirical  model 
ANOVA  tests,  and  the  simplex  optimization  method  successfully  determined  an  optimized 
antenna  design  for  the  microstrip  antenna  configuration  of  Figure  4.1  for  average  antenna 
gain  in  the  angular  regions  in  Equations  (4.1)  and  (4.2).  Validated  by  HFEM  code  compar¬ 
isons  the  optimized  design  methodology  showed  its  strength  in  the  time  savings  observed 
while  performing  the  optimization  procedure.  In  addition,  using  the  cubic  spline  empirical 
model  to  perform  ANOVA  testing  for  the  full  variable  case  quickly  demonstrated  the  an¬ 
tenna  design  sensitivity  to  variable  changes  and  indicated  the  design  areas  requiring  close 
manufacturing  control.  The  next  step  towards  completely  establishing  this  optimized  de¬ 
sign  methodology  for  this  antenna  configuration  is  its  application  to  a  realizable  antenna 
with  accessible  dielectric  values  followed  by  antenna  range  validation. 
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V.  Optimized  Design  Methodology  Implementation 

This  chapter  implements  the  optimized  design  methodology  developed  in  Chapter  IV 
and  uses  a  more  accessible  dielectric  substrate  with  a  relative  permittivity  value  of  4.5. 
Changing  the  substrate’s  dielectric  constant  requires  the  scaling  of  Figure  4.1  which  allows 
for  minor  design  variable  adjustments.  This  implementation  of  the  design  methodology 
uses  a  full  factorial  designed  experiment,  cubic  spline  empirical  modeling,  and  simplex  op¬ 
timization  as  in  Chapter  IV,  but  it  applies  a  combined  optimization  function  that  requires 
minimum  average  antenna  gain  and  50  Ohm  input  impedance.  The  optimization  results 
produce  a  realizable  antenna  design  that  corresponds  well  with  Hybrid  Finite  Element 
Method  (HFEM)  results  and  with  antenna  range  measurements. 

5.1  Experiment  Development 

Since  changing  the  substrate  relative  permittivity  of  Figure  4.1,  scaling  using  Equa¬ 
tion  (5.1)  produced  the  antenna  of  Figure  5.1  with  dimensions  Usted  in  Table  5.1. 

ScaleF  actor  =  (5.1) 

Y  ^new 

In  Equation  (5.1)  Coid  represents  the  substrate  permittivity  of  Figure  4.1,  and  Cnew  rep¬ 
resents  the  substrate  permittivity  of  Figure  5.1.  The  new  antenna  configuration  and  the 
intention  of  optimizing  the  antenna  design  well  within  the  variable  ranges  rather  than  at 
variable  limits  led  to  adjustments  of  the  input  variable  levels  used  in  the  design  method¬ 
ology  development  of  Chapter  IV. 

5AA  Design  Methodology  Input  Variable  Levels.  As  in  Chapter  IV,  the  design 
methodology  input  variables  include: 

1.  Antenna  operating  frequency, 

2.  First  dielectric  superstrate  layer  relative  permittivity, 

3.  Second  dielectric  superstrate  layer  relative  permittivity, 

4.  First  dielectric  superstrate  layer  thickness. 
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1.343  cm 


0.894  cm 


A 

V 


1 1.078  cm - ^ 

6.157  cm  ^ 


Figure  5.1  Scaled  Antenna  Design 


Scale  Factor 

0.94281 

Cavity  Dimensions: 
Length 

Width 

Height 

11.078  cm 

11.078  cm 

2.237  cm 

Patch  Dimensions: 
Length 

Width 

6.157  cm 
6.157  cm 

Patch  Location 

Centered  0.670  cm 
Above  Cavity  Bottom 

Feed  Pin  Location 

0.615  cm  Off  Center 

Substrate  Permittivity 

4.5 

Substrate  Thickness 

0.894  cm 

Table  5.1  Scaled  Constant  Antenna  Design  Parameters 
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5. 1.1.1  Antenna  Operating  Frequency.  As  in  Section  4. 1.3.1,  two  pre¬ 
liminary  HFEM  tests  that  determine  antenna  input  impedances  for  the  scaled  antenna 
configuration  of  Figure  5.1  located  the  frequency  range  of  interest.  Using  first  and  sec¬ 
ond  dielectric  layer  relative  permittivities  of  1.0  each  and  a  characteristic  impedance  of 
50  ohms,  Figure  5.2  shows  the  antenna  input  admittance  response  using  frequencies  from 
1.1  to  1.3  GHz.  Likewise,  using  first  and  second  dielectric  layer  relative  permittivities  of 
10.0  and  a  characteristic  impedance  of  50  ohms.  Figure  5.3  shows  the  antenna  admittance 
response  using  frequencies  from  1.2  to  1.4  GHz.  Since  the  optimized  antenna  design  found 
in  Section  4.4  converged  to  the  upper  frequency  limit  used  in  Chapter  IV,  new  frequency 
levels  are  needed.  Prom  Figures  5.2  and  5.3,  the  new  frequency  levels  include  1,2,  1.25, 
1.3,  1.35,  and  1.4  GHz.  Using  Equation  (5.2), 

/=^.  (5.2) 

where  C  is  the  speed  of  light  and  L  is  the  antenna  cavity  width,  the  new  frequency  level 
range  covers  the  unfilled  antenna  cavity  half  wavelength  resonant  frequency,  1.353  GHz. 

5. 1.1. 2  First  and  Second  Dielectric  Superstrate  Layer  Permittivities.  As 
in  Section  4.1. 3. 2,  using  practical  antenna  design  experience  and  desiring  easily  accessi¬ 
ble  dielectric  materials  led  to  the  selection  of  new  first  and  second  dielectric  superstrate 
layer  relative  permittivities.  In  addition  since  the  optimized  antenna  design  of  Section  4.4 
converged  to  an  upper  limit  for  the  second  dielectric  layer  relative  permittivity,  the  new 
relative  permittivity  values  include  2.0,  4.0,  6.0,  8.0,  and  10.0. 

5. 1.1.3  First  Dielectric  Superstrate  Layer  Thickness.  As  in  Section  4. 1.3.3, 
the  partitioning  property  of  the  HFEM  code  used  in  the  full  factorial  designed  experiment 
limits  the  levels  for  the  first  dielectric  superstrate  layer  to  include  0.2237,  0.4474,  0.6711, 
0.8948,  and  1.1185  centimeters.  Table  5.2  lists  the  new  designed  experiment  variables  and 
their  corresponding  variable  levels. 

5.1.2  Measured  Responses,  The  measured  responses  from  the  full  factorial  de¬ 
signed  experiment  include  both  the  average  antenna  gain  from  Equation  (4.3)  over  the 
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Figure  5.2  Scaled  Antenna  Smith  Chart  Analysis  of  Input  Impedance  Values  With  First 
and  Second  Superstrate  Dielectric  Relative  Permittivities  Equal  to  1.0  (e^i  = 
er2  1.0) 


Frequency  (GHz) 

1.20,  1.25,  1.30 
1.35,  and  1.40 

First  Layer  Relative 
Permittivity 

2,  4,  6, 
8,  and  10 

Second  Layer  Relative 
Permittivity 

2,  4,  6, 
8,  and  10 

First  Layer 

Thickness  (cm) 

0.2237,  0.4474, 
0.6711,  0.8948, 
and  1.1185 

Table  5.2  Scaled  Variable  Antenna  Design  Parameters 


5-4 


Figure  5.3  Scaled  Antenna  Smith  Chart  Analysis  of  Input  Impedance  Values  With  First 
Eind  Second  Superstrate  Dielectric  Relative  Permittivities  Equal  to  10.0  (e^i  = 
er2  =  10.0) 
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Test 

Layer  1 
(cm) 

6rl 

^r2 

..  _ 

Freq 

(GHz) 

Gain 

(dB) 

1 

0.2237 

2 

2 

1.20 

yi 

2 

0.2237 

2 

2 

1.25 

V2 

3 

0.2237 

2 

2 

1.30 

yz 

4 

0.2237 

2 

2 

1.35 

yA 

5 

0.2237 

2 

2 

1.40 

yz 

6 

0.2237 

2 

4 

1.20 

ye 

625 

1.2185 

10 

10 

1.40 

J/625 

Table  5.3  Scaled  Antenna  Pull  Factorial  Design 


Frequency  (GHz) 

1.26,  1.32,  1.38 

First  Layer  Relative 
Permittivity 

3,  7,9 

Second  Layer  Relative 
Permittivity 

3,  7,  9, 

First  Layer 

Thickness  (cm) 

0.2237,  0.6711,  1.1185 

Table  5.4  Scaled  Antenna  Model  Adequacy  Test  Variable  Levels 

angular  region  described  by  Equations  (4.1)  and  (4.2)  and  the  antenna  input  impedance. 
Later,  the  simplex  optimization  method  combines  the  average  antenna  gain  and  the  an¬ 
tenna  input  impedance  to  form  a  combined  optimization  function  to  locate  an  antenna 
design  with  minimum  average  gain  and  an  input  impedance  of  50  ohms. 

5,1,3  Experimental  Design,  Using  a  similar  full  factorial  designed  experiment 
as  Section  4.1.5  with  new  input  variable  levels,  Table  5.3  lists  the  new  four-variable,  five- 
level  designed  experiment.  This  full  factorial  experiment  includes  625  (5^)  system  response 
measurements  for  both  the  average  antenna  gain  and  the  antenna  input  impedance. 


5,2  Average  Antenna  Gain  and  Input  Impedance  Modeling 

Using  the  cubic  spline  multi-dimensional  empirical  modeling  algorithm  introduced  in 
Section  4.2.1,  the  variable  levels  of  Table  5.4,  and  the  cases  listed  in  Table  5.5,  Section  5.2.1 
and  Section  5.2.2  demonstrate  empirical  modeling  adequacies  for  average  antenna  gain  and 
the  antenna  input  impedance  respectively. 
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Test 

Layer  1 
(cm) 

€rl 

^r2 

Preq 

(GHz) 

Gain 

(dB) 

1 

0.2237 

3 

3 

1.26 

yi 

2 

0.2237 

3 

3 

1.32 

y2 

3 

0.2237 

3 

3 

1.38 

2/3 

4 

0.2237 

3 

7 

1.26 

yA 

81 

1.1185 

9 

9 

1.38 

2/81 

Table  5.5  Scaled  Antenna  Model  Adequacy  Test  Cases 

5, 2 A  Average  Antenna  Gain  Model  Adequacy.  Following  the  example  of  Sec¬ 
tion  4.2.3,  this  sections  applies  the  empirical  model  adequacy  tests  of  Section  3.3.4  to 
verify  the  accuracy  of  the  average  antenna  gain  cubic  spline  empirical  model  for  the  new 
scaled  antenna.  Prom  Figure  5.4,  the  empirical  model  varies  over  a  net  range  of  approxi¬ 
mately  2.2  dB,  1.0  dB  greater  than  the  cubic  spline  empirical  model  used  in  Section  4,2.3. 
The  residuals  in  the  normal  probability  plot  of  Figure  5.4  form  less  of  a  straight  line  than 
the  residuals  of  Figure  4.4  and  indicate  less  of  a  good  model  fit.  Fortunately,  the  residual 
mean  and  the  residual  standard  deviation  values  listed  in  Table  5.6  for  the  new  average 
antenna  gain  empirical  model  closely  correspond  to  those  of  Table  4.6.  Recall  that  the  new 
scaled  antenna  employs  higher  dielectric  permittivities  than  the  previous  antenna  config¬ 
uration.  Therefore,  the  increased  variability  of  the  empirical  model  most  likely  originates 
from  the  larger  superstrate  and  substrate  permittivities.  Figure  5.5  indicates  somewhat  of 
a  normal  random  distribution  of  model  residual  error  with  very  close  average  gain  pred¬ 
ications  between  -2.0  dB  and  2.0  dB.  Figures  5.6,  5.7,  5.8,  and  5.9  show  the  effects  of 
variable  levels  on  the  average  antenna  gain.  Prom  Figure  5.6,  the  test  case  frequencies, 
1.32  and  1.38  GHz,  strongly  affected  the  average  gain  since  both  reside  near  the  unfilled 
cavity  half  wavelength  resonant  frequency,  1.353  GHz,  Furthermore,  from  Figure  5.7  high 
relative  permittivity  values  from  the  first  dielectric  layer  caused  increased  model  residual 
values,  while  lower  relative  permittivity  values  from  the  second  dielectric  layer  seen  in 
Figure  5.8  caused  increased  model  residual  values.  Worth  noting  here,  the  thickness  of  the 
second  dielectric  superstrate  layer  directly  depends  on  the  thickness  of  the  first  dielectric 
superstrate  layer.  Prom  Figure  5.1,  when  the  first  superstrate  layer  thickness  equals  0.6711 
cm,  the  second  superstrate  layer  thickness  also  equals  0.6711  cm.  Figure  5.9  indicates  that 
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Normal  Probability  Plot 


Average  Antenna  Gain  Residuals  (dB) 

Figure  5.4  Average  Antenna  Gain  Empirical  Model  Normal  Probability  Plots  of  Resid¬ 
uals 

Residual  Mean  -0.1431  dB 

Residual  Standard  Deviation  0.3567  dB 

Table  5.6  Average  Antenna  Gain  Empirical  Model  Residual  Mean  And  Residual  Stan¬ 
dard  Deviation  Values 

when  the  first  superstrate  layer  thickness  equals  0.6711  cm,  the  empirical  model  remains 
accurate  with  low  model  residuals.  Despite  the  increased  variability  of  the  average  an¬ 
tenna  gain  empirical  model,  the  model  still  successfully  predicts  antenna  gain  trends  and 
still  provides  a  useful  gain  prediction  tool. 

5.2.2  Antenna  Input  Impedance  Model  Adequacy.  As  in  Section  5.2.1,  this  sec¬ 
tions  employs  the  empirical  model  adequacy  tests  of  Section  3.3.4  to  verify  the  accuracy 
of  the  antenna  input  impedance  cubic  spline  empirical  model  for  both  real  and  imaginary 
impedances.  Prom  Figure  5.10,  neither  the  real  nor  the  imaginary  impedance  empirical 
models  precisely  predict  antenna  input  impedance.  Worth  noting,  the  real  and  imagi¬ 
nary  impedance  residual  means  and  residual  standard  deviations  listed  in  Table  5.7  closely 
match  those  listed  in  Table  4.14  for  the  previous  configuration  but  still  remain  large.  Fig¬ 
ures  5.11,  5.12,  5.13,  5.14,  and  5.15  indicate  the  wide  variability  of  the  antenna  input 
impedance  models  and  show  no  sign  of  strong  antenna  impedance  influences  do  to  a  single 


Cubic  Spline  Fitted  Values  (dB) 

Figure  5.5  Average  Antenna  Gain  Empirical  Model  Residuals  Versus  Fitted  Values 


Antenna  Frequency  (GHz) 


Figure  5.6  Average  Antenna  Gain  Empirical  Model  Residuals  Versus  Frequency  Values 
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Figure  5.7  Average  Antenna  Gain  Empirical  Model  Residuals  Versus  First  Layer  Rela¬ 
tive  Permittivity  Values 


Figure  5.8  Average  Antenna  Gain  Empirical  Model  Residuals  Versus  Second  Layer  Rel¬ 
ative  Permittivity  Values 
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Figure  5.9  Average  Antenna  Gain  Empirical  Model  Residuals  Versus  First  Layer  Thick¬ 
ness  Values 


Mean  Real  Impedance  Residual 

Real  Impedance  Residual  Standard  Deviation 

1.3091  ohms 
4.1565  ohms 

Mean  Imaginary  Impedance  Residual 

Imaginary  Impedance  Residual  Standard  Deviation 

1.5775  ohms 
6.5177  ohms 

Table  5.7  Real  and  Imaginary  Antenna  Input  Impedance  Empirical  Model  Mean  Resid¬ 
ual  and  Residual  Standard  Deviations  Values 

variable  or  variable  level.  Fortunately,  like  the  average  antenna  gain  empirical  model  of 
Section  5.2.1  the  impedance  empirical  models  still  predict  impedance  trends  and  provide 
a  useful  tool  for  antenna  design  optimization. 
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Probability  Probability 
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Figure  5.10  Antenna  Input  Impedance  Empirical  Model  Normal  Probability  Plots  of 
Residuals 


5-12 


Residuals  (ohms) 


10  15  20  25  30 

Real  Impedance  Fitted  Values  (ohms) 


Imaginary  Impedance  Fitted  Values  (ohms) 


Figure  5.11  Antenna  Input  Impedance  Empirical  Model  Residuals  Versus  Fitted  Values 
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Figure  5.13  Antenna  Input  Impedance  Empirical  Model  Residuals  Versus  First  Layer 
Relative  Permittivity  Values 
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Figure  5.14  Antenna  Input  Impedance  Empirical  Model  Residuals  Versus  Second  Layer 
Relative  Permittivity  Values 
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Imaginary  Residuals  (ohms) 


Siy 


Frequency 

1.30  GHz 

First  Layer  Relative  Permittivity 

6.0 

Second  Layer  Relative  Permittivity 

10.0 

First  Layer  Thickness 

0.4474  cm 

Average  Gain  (HFEM) 

-4.7748  dB 

Input  Impedance  (HFEM) 

19.5875  -1-  J67.6017  ohms 

Optimized  Solution 

25.638 

Table  5.8  Optimization  Starting  Point 


5.3  Average  Antenna  Gain  and  Antenna  Input  Impedance  Optimization 

Since  the  goal  of  this  chapter  is  the  design  of  a  realizable  antenna  for  the  configuration 
in  Figure  5.1,  the  optimization  requirements  include: 

1.  Minimum  average  antenna  gain  in  the  angular  region  described  in  Equations  (4.1) 
and  (4.2). 

2.  Input  impedance  of  50  ohms. 

Using  the  simplex  optimization  method  as  in  Section  4.4,  the  empirical  models  of  Sec¬ 
tions  5.2.1  and  5.2.2,  and  the  combined  optimization  function  of  Equation  (5.3), 

solution  ~  averagegain  +  a6s(50  —  real{impedance))  (5.3) 

the  optimization  procedure  searches  through  antenna  radiation  patterns  like  those  of  Fig¬ 
ures  5.16  and  5.17  to  satisfy  the  optimization  requirements.  Using  the  optimization  starting 
point  of  Table  5.8  with  its  radiation  pattern  shown  in  Figure  5.18,  the  simplex  optimization 
method  converged  to  the  optimized  antenna  design  location  listed  in  Table  5.9  with  huge 
improvements  over  the  starting  point  design.  Figure  5.19  provides  the  optimized  design 
radiation  pattern,  while  Figures  5.20,  5.21,  5.22,  5.23,  5.24,  and  5.25  show  the  combined 
optimization  response  as  variable  pairs  approach  their  optimized  values. 
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2.5 


Linear  "X"  Directed  Gain  Linear  "Y"  Directed  Gain 


Figure  5.16  Antenna  Radiation  Pattern  (Frequency  =  1.2  GHz,  First  Superstrate  Di¬ 
electric  Layer  Permittivity  =:  Second  Superstrate  Dielectric  Layer  Relative 
Permittivity  —  2.0  ) 


Figure  5,17  Antenna  Radiation  Pattern  (Frequency  —1.4  GHz,  First  Superstrate  Di¬ 
electric  Layer  Permittivity  =  Second  Superstrate  Dielectric  Layer  Relative 
Permittivity  =  10.0  ) 
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Frequency 

1.353  GHz 

First  Layer  Relative  Permittivity 

6.010 

Second  Layer  Relative  Permittivity 

9.966 

First  Layer  Thickness 

0.673  cm 

Elapsed  Optimization  Time 

2  mins 

Average  Gain  (Empirical  Model) 

-6.501  dB 

Input  Impedance  (Empirical  Model) 

50.00  -|-  j 38.086  ohms 

Average  Gain  (HFEM) 

-6.794  dB 

Input  Impedance  (HFEM) 

48.611  +  j36.232  ohms 

Gain  Accuracy 

4.3  percent 

Impedance  Accuracy 

3.7  percent  (magnitude) 

Weighted  Optimization  Value 

-6.501  (Empirical  values) 

Weighted  Optimization  Value 

-5.405  (HFEM  values) 

Table  5.9  Optimized  Antenna  Design 
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Combined  Optimization  Values  Linear  "Z"  Directed  Gain 


I  20s  . 
E 

■S.  10x..- 
o 


Second  layer  permittivity 


Frequency  (GHz) 


Figure  5.21  Second  Layer  Relative  Permittivity  and  Frequency  Effects  Near  the  Opti¬ 
mized  Location 


Frequency  (GHz) 


Figure  5.22  First  Layer  Thickness  and  Frequency  Effects  Near  the  Optimized  Location 
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40 


10  5 


Second  layer  permittivity  First  layer  permittivity 

Figure  5.23  First  and  Second  Layer  Relative  Permittivity  Effects  Near  the  Optimized 
Location 


Figure  5,24  First  Layer  Thickness  and  Second  Layer  Relative  Permittivity  Effects  Near 
the  Optimized  Location 
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Figure  5.25  First  Layer  Relative  Permittivity  and  Thickness  Effects  Near  the  Optimized 
Location 
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Scale  Factor 

3.5236 

Frequency 

4.768  GHz 

First  Superstrate  Layer  Relative  Permittivity 

,6.15 

First  Superstrate  Layer  Thickness 

0.075  inches  (0.1905  cm) 

Second  Superstrate  Layer  Relative  Permittivity 

10.2 

Second  Superstrate  Layer  Thickness 

0.075  inches  (0.1905  cm) 

Substrate  Relative  Permittivity 

4.5 

Substrate  Thickness 

0.10  inches  (0.254  cm) 

Cavity  Dimensions: 

Length 

1.2375  inches  (3.1433  cm) 

Width 

1.2375  inches  (3.1433  cm) 

Height 

0.250  inches  (0.635  cm) 

Patch  Dimensions: 

Length 

0.550  inches  (1.397  cm) 

Width 

0.550  inches  (1.397  cm) 

Patch  Location 

Centered  0.075  inches  (0.1905  cm) 

Above  Cavity  Bottom 

Feed  Pin  Location 

0.0687  inches  (0.175  cm)  Off  Center 

Predicted  average  gain  (HFEM) 

-6.582  dB 

Predicted  input  impedance  (HFEM) 

45.324  -1-  y25.694  ohms 

Measinred  Input  Impedance 

51.439  4-i3.045  ohms 

Table  5.10  Constructed  Antenna  Design  Parameters 


5,4  Compact  Range  Validation  of  Designed  Cavity-Backed  Microstrip  Antenna  With  Di¬ 
electric  Overlays 

Using  the  optimized  design  of  Table  5,9,  the  completed  optimized  design  methodology 
requires  a  constructed  antenna  for  methodology  verification.  Unfortimately,  exact  dielec- 
trie  materials  with  the  permittivity  values  of  Table  5.9  are  not  readily  available.  Therefore, 
the  available  dielectric  materials  for  superstrate  layers  one  and  two  are  6.15  and  10.2  re¬ 
spectively.  Table  5.10  lists  the  antenna  design  parameters  including  the  new  antenna  scale 
factor,  the  new  antenna  dimensions,  and  predicted  antenna  responses.  Figure  5.26  shows 
the  new  scaled  antenna  design  and  Figure  5.27  shows  the  HFEM  predicted  radiation  pat¬ 
tern  with  a  predicted  average  antenna  gain  of  -6.582  dB  and  a  predicted  input  impedance 
of  45.324  +  j25.694  ohms.  The  resulting  differences  between  the  predicted  average  antenna 
gain  and  input  impedance  values  listed  in  Table  5,10  for  the  constructed  antenna  and  Ta¬ 
ble  5.9  for  the  optimized  antenna  design  are  directly  due  to  using  the  available  dielectric 
materials  for  the  first  and  second  superstrate  layers  rather  than  the  exact  optimized  dielec- 
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Figure  5.26  Constructed  Antenna 


trie  values.  In  fact,  the  measured  input  impedance  from  the  constructed  antenna  measured 
51.439  +j3.045  ohms,  where  the  difference  is  due  to  the  antenna  radiation  aperture’s  feed 
pin  location. 

Measurement  of  the  constructed  antenna  using  a  compact  anechoic  antenna  chamber 
to  verify  the  optimized  design  methodology  and  HFEM  predictions  produced  the  electric 
field  patterns  shown  in  Figures  5.28,  5.29,  5.30,  and  5.31  where  the  smooth  curves  are 
predicted  electric  field  patterns  and  the  rough  curves  are  the  measured  electric  field  pat¬ 
terns.  In  general,  the  measured  fields  closely  match  the  predicted  fields  with  only  minor 
differences.  In  Figures  5.29  and  5.31,  two  differences  exist.  The  first,  a  small  hump  in 
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Figure  5.27  Predicted  Radiation  Pattern  for  the  Constructed  Antenna  Design 

the  electric  field  pattern  near  the  zero  degree  angle,  is  caused  by  uneven  mounting  of  the 
constructed  antenna  into  a  measurement  test  body.  The  second,  electric  field  round  off  at 
the  field  extremities,  is  caused  by  the  finite  ground  plane  used  for  antenna  measurements 
instead  of  the  ideal  infinite  ground  plane  used  in  antenna  field  predictions.  In  Figure  5.30, 
the  expected  electric  field  was  approximately  -80  dB,  where  the  measured  electric  field  was 
in  the  general  range  between  -20  dB  and  -40  dB.  The  difference  is  attributed  to  cross  po¬ 
larization  effects  due  to  misalignment  of  the  test  body  during  field  pattern  measurements. 


5.5  Closing  Comments 

Using  the  optimized  design  methodology  of  Chapter  IV  and  a  combined  optimization 
constraint  equation  composed  of  minimum  average  antenna  gain  and  a  target  input  im¬ 
pedance  of  50  ohms,  this  chapter  located,  built,  and  verified  a  realizable  optimized  antenna 
design.  Despite  increased  variations  in  the  average  antenna  gain  and  input  impedance 
empirical  models,  the  optimized  antenna  design  compared  well  with  HFEM  results  and 
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measured  results.  The  constructed  antenna  used  to  verify  the  optimized  antenna  design 
used  existing  dielectric  materials  that  closely  matched  the  optimized  dielectric  materials 
and  produced  electric  field  patterns  that  closely  resembled  the  predicted  electric  fields.  As 
a  result  this  RSM  based  design  methodology  is  an  excellent  design  tool  for  this  antenna 
configuration. 
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VL  Conclusion 


This  thesis  introduced  a  relatively  new  approach  to  the  design  of  cavity-backed  mi¬ 
crostrip  antennas  that  uses  RSM  techniques  to  determine  an  optimized  antenna  design. 
Section  6.1  states  some  specific  conclusions  found  using  the  design  methodology  of  Chap¬ 
ter  IV  and  its  implementation  in  Chapter  V  and  summarizes  some  of  the  key  impacts  found. 
Section  6.2  lists  some  areas  where  further  use  or  extension  of  this  design  methodology  is 
possible. 

6.1  Summary  and  Conclusions 

Prom  this  research,  several  key  conclusions  are  drawn  that  highlight  the  utility  of 
this  design  methodology  and  the  individual  RSM  techniques  used  to  design  an  optimized 
cavity-backed  microstrip  antenna  with  dielectric  overlays.  The  conclusions  include: 

1.  The  RSM  based  design  methodology  quickly  and  accurately  predicted  an  optimized 
antenna  design  for  this  particular  antenna  configuration  where  the  design  method¬ 
ology  results  closely  agreed  with  HFEM  comparison  results  and  antenna  range  mea¬ 
surements.  In  general,  all  optimized  results  were  within  4.0  percent  of  the  HFEM 
results,  and  the  predicted  electric  field  patterns  closely  resembled  actual  antenna 
range  field  pattern  measurements  with  only  slight  range  measurement  differences. 

2.  The  full  factorial  designed  experiment  along  with  the  design  variables  and  variable 
levels  suited  this  design  methodology  well.  The  full  factorial  designed  experiment 
allowed  for  the  drastically  changing  average  antenna  gain  and  antenna  input  im¬ 
pedance  values  to  be  sampled  throughout  the  entire  region  of  design  interest  where 
each  design  variable  strongly  impacted  the  overall  design  problem.  Due  to  the  cav¬ 
ity  configuration  of  this  antenna,  antenna  operating  frequency  directly  infiuenced 
the  resulting  antenna  performance,  and  due  to  dielectric  loading  effects,  the  rela¬ 
tive  permittivity  variables  of  superstrate  layers  one  and  two  as  well  as  the  thickness 
variable  of  superstrate  layer  one  directly  impacted  the  resulting  antenna  gain,  input 
impedance,  and  antenna  radiation  pattern  shape. 
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3.  The  cubic  spline  empirical  model  used  throughout  this  research  quickly  and  accu¬ 
rately  determined  average  antenna  gain  with  relatively  small  residual  means  and 
residual  standard  deviations  but  had  difficulty  accurately  modeling  and  predicting 
antenna  input  impedance.  As  mentioned  in  Chapter  II,  dielectric  loading  reduces 
the  antenna  input  impedance  bandwidth  and  causes  the  impedance  bandwidth  to 
become  even  more  restrictive  than  the  frequency  bandwidth  [17,18].  In  other  words, 
slight  changes  in  the  relative  permittivity  of  either  superstrate  layer  one  or  layer 
two  greatly  impacted  the  resulting  antenna  input  impedance  for  this  antenna  design 
causing  the  input  impedance  to  have  sharp  changes  in  value. 

4.  The  ANOVA  tests  performed  in  Chapter  IV  not  only  indicated  the  impacts  that 
each  design  variable  and  each  variable  interaction  factor  had  on  the  resulting  average 
antenna  gain  performance  for  the  antenna  of  Figure  4.1  but  also  provided  a  design  tool 
indicating  where  design  sensitivities  were  the  strongest.  Using  ANOVA  testing  in  this 
manner  shows  manufacturers  and  engineers  where  manufacturing  tolerances  should 
be  the  tightest.  Prom  Section  4.3.2,  the  antenna  operating  frequency  and  the  second 
dielectric  superstrate  layer  relative  permittivity  factors  were  the  two  design  factors 
with  the  largest  impact  on  the  average  antenna  gain  therefore  requiring  the  tightest 
manufacturing  control.  This  type  of  statistical  analysis  is  only  feasible  because  of 
the  efficiency  of  the  cubic  spline  empirical  model. 

5.  The  application  of  the  simplex  optimization  method  to  the  cubic  spline  empirical 
models  in  both  Chapter  IV  and  Chapter  V  quickly  determined  an  optimized  antenna 
design  that  compared  well  with  verification  results.  In  Chapter  IV  during  the  design 
methodology  development  phase,  the  simplex  optimization  method  used  the  cubic 
spline  empirical  model  to  determine  an  optimized  antenna  design  in  an  elapsed  time 
interval  of  approximately  one  minute  using  180  design  iterations.  The  HFEM  opti¬ 
mized  design  results  verification  test  took  approximately  three  hours.  Assuming  each 
simplex  iteration  would  have  taken  as  long,  it  would  have  taken  23  days  to  come  up 
with  the  optimized  design  using  HFEM  testing  rather  than  empirical  modeling.  The 
time  savings  found  using  empirical  model  based  optimization  is  the  power  of  RSM 
and  designed  experiments. 
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6, 2  Recommendations 


The  results  found  in  this  research  indicate  that  the  design  methodology  developed 
here  is  extendible  to  other  electromagnetic  design  problems  with  the  proper  selection  of 
design  variables,  designed  experiment,  and  empirical  model.  To  build  confidence  in  the 
validity  of  using  this  design  methodology,  several  areas  for  continued  research  and  design 
are  listed  below: 

1.  Optimize  for  a  maximum  antenna  gain  using  the  antenna  configuration  of  Figure  4.1 
in  a  designated  angular  region  for  a  more  directive  antenna  radiation  pattern. 

2.  Explore  the  usefulness  of  other  less  time-expensive  designed  experiments  such  as  the 
central  composite  design  or  the  Latin  squares  design  since  the  full  factorial  designed 
experiment  required  such  a  large  number  of  antenna  response  measurements. 

3.  Compare  this  design  methodology  with  other  design  approaches  in  recent  antenna 
design  literature  such  as  those  using  genetic  algorithms  or  quadratic  optimization. 

4.  Perform  rigorous  ANOVA  test  comparisons  between  the  HFEM  and  the  empirical 
model  results  for  design  parameters  other  than  frequency  to  build  confidence  in  the 
validity  of  using  an  empirical  model  in  ANOVA  testing. 

5.  Apply  the  design  methodology  to  other  antenna  designs  such  as  a  spiral  antenna 
design. 

6.  Apply  the  design  methodology  to  other  electromagnetic  problems  such  as  periodic 
surfaces. 
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